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1.0 INTRODUCTION 


Over the course of the Magsat project there has been much debate 
about the importance of large-scale remanent magnetization in the crust of the 
continents, and speculation about whether Magsat data can detect the phenom- 
enon if it exists. These questions were the subject of the present investiga- 
tion. 


Magnetic anomaly maps derived from Pogo data (Regan et al , 1975; 
Mayhew, 1982) and Magsat data (Langel et al , 1982; Mayhew and Galliher, 1982) 
show a prominent anomaly over Kentucky and Tennessee (Figure 1). Equivalent 
layer magnetization models derived by inversion of such data (e.g. Mayhew and 
Galliher, 1982) indicate an extremely magnetic source region centered in 
Kentucky (Figure 2). The size of the magnetization anomaly suggests that this 
is the most important large-scale concentration of magnetic material in the 
crust of the continental U.S., yet there is no obvious direct expression of 
the source at the surface. 

Mayhew et al (1982) found that a prominent elongate gravity anomaly 
occurs at the center of the source region indicated by magnetization models 
based on satellite data. A long wavelength aeromagnetic anomaly is directly 
associated with the gravity anomaly, although its form is largely masked by 
local near-surface magnetic anomalies. Using limited crustal seismic refrac- 
tion data for constraint, the above authors produced a simple model, repre- 
sented as two dimensional cross-sections, which accounts for the gravity and 
associated aeromagnetic anomalies. The model is in the form of an elongate 
body which is anomalously dense and magnetic extending through most of the 
thickness of the crust. On the basis of several lines of evidence, it was 
interpreted as a large mass of mafic intrusive rock of late Precambrian age, 
and was termed the "Kentucky body". Keller et al (1976) considered the 
gravity high to be part of a more extensive belt which they termed the "East 
Continent Gravity High", and interpreted as the expression of a late 
Precambrian rift zone. 
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For the first part of this study a three-dimensional model of the 
Kentucky body was developed to fit surface gravity and long wavelength aero- 
magnetic data. Magnetization and density parameters for the model are much 
like those of Mayhew et al (1982). The magnetic anomaly due to the model at 
satellite altitude is shown to be much too small by itself to account for the 
anomaly measured by Magsat. It is demonstrated that the source region for the 
satellite anomaly is considerably more extensive than the Kentucky body sensu 
stricto . The extended source region is modeled in the second part of the 
study, first using prismatic model sources and second using dipole array 
sources. Magnetization directions for the source region found by inversion of 
various combinations of scalar and vector data are found to be close to the 
main field direction, implying the lack of a strong remanent component. It is 
shown by simulation that in a case (such as this) where the geometry of the 
source is known, if a strong remanent component is present its direction is 
readily detectable, but by scalar data as readily as vector data. 
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2.0 LOCAL MODEL OF KENTUCKY BODY 


Figure 3a is a local Bouguer gravity anomaly map for the area of the 
Kentucky body, while Figure 3b is the aeromagnetic map for the corresponding 
area. Figure 4 is the outline of a vertical -sided prismatic model body 
constructed to produce gravity and magnetic anomalies giving a gross fit to 
those observed. The model is divided into three parts. The top of the main 
(central) part is placed at 6 km below sea level; constraint for this is the 
refraction line of Warren (1968). The bottom is poorly constrained, but is 
placed at 42 km depth, about 8 km beneath the local Moho, giving a negative 
density contrast for this section which may account for the negative side- 
lobes in the computed gravity (Figure 5). The central body is approximately 
in isostatic equilibrium. The refraction line of Borcherdt and Roller (1966) 
extending across the south end of the gravity high indicates an abrupt drop in 
depth to the anomalous body from 6 km in the central part to 16 km in the 
southern part. This was used as a constraint on depth to top for the southern 
part of the model body. Calculations suggest that the southern part does not 
extend as deep below the local Moho as the central part, but there is no good 
constraint on this. For the final model, bottom was taken to be at Moho 
depth. The form of the north end of the gravity anomaly suggests a deepening 
of the top of the source. The Irvine-Paint Creek fault zone bounds the 
gravity anomaly on the north. Faulting is down-dropped to the north, consis- 
tent with the above inference. It is assumed that depth to the top of the 
body is 16 km, as in the south; this is consistent with the cross-section of 
Ammerman and Keller (1979) just to the east. Bottom is taken to be at Moho 
depth. The margins of the model body were placed by trial and error, using 
the gravity gradients as a guide, until a reasonably good fit was obtained. 

No attempt was made to fit the surrounding anomalies, for which there are no 
constraints. While there is not enough information to constrain the details 
of the geometry of the source, the gross geometry and density are well 
determined. 

In modeling the source, the body was first constructed to fit the 
gravity anomaly (Figure 5) using Plouff's (1976) algorithm. Then the assumed 
magnetization for the body was found such that the computed magnetic anomaly 
amplitude was in agreement with the long wavelength part of that observed 
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(Figure 6). Again, it is not possible to fit the detailed local anomalies, 
but the mean magnetization of the model body itself (5.2 A/m) is well deter- 
mined, certainly within 10%. Magnetic model computations also used Plouff's 
approach. 


Finally, the magnetic anomaly due to the Kentucky body at satellite 
altitude (325 km) was computed, using the inferred magnetization value (Figure 
7). Clearly, the computed anomaly is too small by a factor of about three to 
account for the observed anomaly (Figure 1). 
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3.0 REr7l0NAL MAGNETIZATION MODELS 


The reason why the Kentucky body in itself cannot be the sole source 
of the satellite observed magnetic anomaly can now be easily seen from the 
recently-published aeromagnetic map of the U.S. (Zietz, 1982). This map shows 
clearly that the magnetic source region is considerably more extensive than 
the Kentucky body itself. An indication of the extent of this region is 
indicated in Figure 8. The relation of the Kentucky body itself in relation 
to regional tectonic elements can be seen in Figure 9. The extended source 
region was modeled in two ways, first using prismatic sources, and second 
using an array of dipole sources. 

Examination of the U.S. aeromagnetic map indicates that the eastern 
mid-continent is a regionally magnetic high area relative to the eastern 
seaboard region. The gradient separating the two regions is a long straight 
zone known as the New York-Alabama lineament (King and Zietz, 1978), which 
passes just to the east of the Kentucky body. Expression of this lineament is 
present in the U.S. magnetization map (Figure 2), although it is distorted by 
the presence of the Kentucky source region. 

In our modeling we assumed three simple regions, the first and second 
regions representing those parts to the northwest and southeast, respectively, 
of the New York-Alabama lineament, the third representing the extended 
Kentucky source region itself. The purpose of modeling the first and second 
regions along with the Kentucky region was to remove, to the extent possible, 
the biasing effect due to the difference in level between the first two 
regions. 


3.1 Prismatic Models 


These models used the formulation of Plouff (1976) for the magnetic 
anomaly due to vertical polygonal prisms. Prisms were arbitrarily made 40 km 
thick, i.e. comparable with the thickness of the whole crust, like the 
Kentucky body model described in Section 2.0. Calculations with these models 
involve an inherent flat-earth assumption, but for the limited area considered 
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differences with more rigorous spherical -earth models are very minor. For the 
area treated, three polygonal prismatic model sources were used. Two of the 
prisms are very large and have a common boundary, which is the New York- 
Alabama lineament (Figures 10-12). The third prismatic element is intended to 
model the extended source region of the Kentucky anomaly. This third region 
was modeled in three different ways; these are indicated by the blackened 
areas of Figures 10-12. The first model (Figure 10) represents the most 
restricted geographic source area, corresponding to the highest amplitude 
anomalies seen in the U.S. aeromagnetic map. It is actually two small, 
separate, but nearby, sub-regions which collectively give rise to a single 
observed anomaly. The second model (Figure 11) represents the largest area 
which can reasonably bound the Kentucky source region, again based on the 
U.S. magnetic anomaly map. The third area (Figure 12) is somewhat more 
restricted in area than the second, and probably represents the most realistic 
estimate of the boundaries of the source region, in as much as it avoids the 
linear anomalies on the east which are directly associated with the New York- 
Alabama lineament, rather than the Kentucky source region. 

With the source geometries thus defined, a series of computer runs 
were made in which various combinations of vector component and scalar data 
from an equivalent source reduction of Magsat data (Figure 1) were inverted to 
magnetization solutions for the three regions (two large, one small), and for 
the largest and smallest geometries described above for the Kentucky source 
region. Two types of solutions were obtained. In the first, magnetization 
directions were constrained to be coincident with the main field direction 
("induced" magnetization), and magnitudes only were solved for. In the 
second, magnetization directions were left unconstrained, and solutions were 
found for the magnetization components; this constituted a test for rema- 
nence. Results for a selection of key runs are summarized in Tables 1 and 2. 
These are for 1) input scalar data only, 2) input vector data only, 3) input 
both scalar and vector data. 

Solutions for the two large source regions are not considered partic- 
ularly significant, but it is noted that the second region (southeast) is 
generally more negative than the first (northwest), whether in magnitude or in 
vector orientation, so that incorporation of the two regions in the solutions 
seems to have accounted for some of the regional bias. 
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Solutions for the third (Kentucky source) region are quite consistent 
within each geometry. For the small geometry (Table 1, Figure in) the magnet- 
ization vector magnitude for constrained and unconstrained cases falls between 
7 and 9 A/m, the values being somewhat higher for the unconstrained case. For 
the unconstrained case, the angle between the magnetization vector and the 
main field direction varies from 7 ± 1 degrees to 14 ± 2 degrees, which is not 
significantly different from the main field direction. While the statistics 
are somewhat better where vector data is input, this may simply be the result 
of effectively using more data. 

For the case of the larger, more detailed, prism (Table 2, Figure 
11), angles between the magnetization vector and the main field direction for 
the unconstrained case are quite similar to those for the simple geometry. 
Source volume is larger so we expect magnetization values to be smaller; they 
fall in the range 3-3.5 A/m. 

The anomaly in the total field for the "best" solutions for the 
Kentucky source region of Figure 11 alone for both constrained and uncon- 
strained cases was computed at the altitude of the input data set, 325 km 
(Figures 13 and 14, respectively). The anomalies are dipole-like, with 
slightly different orientations. Note the low on the north side of each 
anomaly. The position of the low is in better agreement with the position of 
a low seen in the same area in the input anomaly data (Figure 1) for the 
constrained case than for the unconstrained case. This perhaps adds some 
support to the idea that magnetization in the source region is by induction in 
the main field. 

A single run was made for the third geometry for the Kentucky source 
region, with magnetization of the source constrained to lie in the main field 
direction. The magnetization solution for this geometry was 4.2 A/m. Note 
that this value is in good agreement with that determined for the Kentucky 
body (5.2 A/m) as described in Section 2.0. This is an interesting result, 
because it suggests that the depth extent for the whole of the extended source 
region is comparable with that of the Kentucky body, i.e. most of the crustal 
thickness. While there are no constraints on this thickness, and 40 km was 
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assumed in the models, the depth extent could not be much less without associ- 
ated magnetization values becoming unreasonably large. The anomaly in the 
total field (AB) due to this geometry for the Kentucky source region alone is 
show in Figure 15. When this is subtracted from the original input aB data 
(Figure 1), an anomaly data set free of the influence of the Kentucky source 
region should result. This is shown in Figure 16. The large "bul 1 ' s-eye" 
(which is the "Kentucky anomaly" of Figure 1) is gone and the remaining 
smaller highs are interpreted to be associated with the New York- Alabama 
1 ineament. 

Finally, as a check on detectability of remanence, we computed the 
magnetic anomaly due to the Kentucky source region alone with an assumed 
direction of magnetization 90° away from the main field direction, added it to 
the residual data set (Figure 16), and then attempted to recover the direction 
using the inversion program. The assumed direction was recovered exactly. 

Our interpretation of the above results corresponds to conventional wisdom: 
where the geometry of the source is reasonably well known, its magnetization 
vector can be found by least squares estimation using one or any combination 
of vector and scalar data types. 

3.2 Dipole Array Models 


These models utilized a new approach to spherical earth equivalent 
dipole modeling of regional sources by constraining all dipoles falling within 
specified regions to adjust together in a least squares solution. The method 
is referred to as mosaic dipole modeling in that a region of interest is 
divided into a set of mosaic subregions in which the dipole arrays are 
constrained. These mosaic subregions are specified by geologic structural 
considerations. A more detailed description of the method is given in 
Appendix A. 


The purpose of utilizing two different approaches for our regional 
magnetization modeling is to provide greater credence and support for conclu- 
sions drawn from the resulting solutions. As with the prismatic models of 
Section 3.1, three mosaic regions were utilized to describe the area treated. 
The dipole grid array and the corresponding mosaic regions are displayed in 
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Figure 17. Regions I and II are large and have a common boundary approximated 
by the New York-Alabama lineament. Region III models the Kentucky body. To 
more closely correspond to the structure represented by the prism of Figure 
10, the three dipoles comprising the Kentucky body were shifted slightly away 
from the grid points displayed in Figure 17. In particular, the western-most 
dipole was shifted .1° to the north and .4° to the east, the northern-most 

dipole was shifted .4° to the south and .5° to the east, and the remaining 

dipole was shifted .4° to the south and .4° to the east. 

With this source geometry, a series of computer runs were made with 

the same combinations of data and solution types as described in Section 3.1. 
The results are presented in Table 3 and show a satisfying consistency with 
the results of the prismatic models. As with the prismatic models, the solu- 
tions for the two large mosaic regions are not considered significant, but do 
show the same more negative structure in the southeast mosaic. The solutions 
for the Kentucky body produce magnetization vector magnitudes very similar to 
those of Section 3.1, while the angle between the magnetization vector and the 
main field direction for the source component solutions varies from 14 ± 5 
degrees to 17 ± 1 degrees. 
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4.0 CONCLUSIONS 


Based on the results of the local model of the Kentucky body and the 
regional magnetization models, we reach the following conclusions: 

1. The Kentucky body alone cannot account for the magnetic anomaly 
measured by Magsat. 

2. The anomaly measured by Magsat can be accounted for by a more 
extensive magnetic source region defined by aeromagnetic data, 
with magnetization (4.2 A/m) comparable with that determined for 
the Kentucky body (5.2 A/m). 

3. The magnetization value for the extended area suggests that 
magnetic material is distributed through most of the crustal 
thickness, as it is for the Kentucky body itself; however, a 
gravity anomaly is not associated with the extended source 
region. 

4. The direction of magnetization determined for the extended source 
region is near the main field direction, suggesting that a strong 
remanent magnetization is lacking. 

5. Where the geometry of the source is known, magnetization 
direction can be found by inverse methods, but equally well with 
scalar data as with vector data. 

6. Results from dipole array modeling agree with results from 
prismatic source modeling; the two very different approaches 
provide a mutual check. 
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APPENDIX A. MOSAIC DIPOLE MODEL 

This technique represents a new method for regional magnetization 
modeling using a constrained dipole equivalent source approach. The method 
consists of using a spherical earth equivalent dipole source model with the 
dipoles within geologically specified mosaic regions constrained to adjust as 
a fixed entity (i.e. all dipoles within the region are constrained to the same 
magnitude and direction by mathematical equations). A least squares 
estimation algorithm best fits the anomaly data while adjusting the 
constrained dipole regions. The program will operate in two modes: 

a) adjust source magnitudes with the directions forced to lie in the main 
field direction 

b) adjust both source magnitudes and direction. 

Data input to the software is any combination of aB^ , ABg , aB^ or aB . 

AAA 

The coordinate system utilized is the spherical r , 0 , (|) system. 

The mathematical description of the least squares algorithm with the 
constrained mosaic regions is given in Appendx A.l, while the mathematical 
derivation of the dipole source function is presented in Appendix A. 2. A 
description of program input is given in Appendix A. 3, while a source listing 
is presented in Appendix A. 4. 
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APPENDIX A.l Least Squares Analysis with Constrained Regions 

Suppose we are given the linear system 

y = Ap + V (Al.l) 

where y is a vector of observations of dimension m , p is a vector of 
parameters of dimension n to be estimated, A is the (m x n) matrix of 
partial derivatives of the modeled observations with respect to the 
parameters, and v is a vector of observation errors of dimension m with 

A 

zero mean, E(v) = D . Then the least squares solution of p of Equation 
(Al.l) is chosen to minimize the square of the observation errors, 

J(p) = (y-AP)"^ (y-AP) = v"^v . (A1.2) 


In the situation where observations of different noise characteristics are 
involved, a solution of the weighted least squared problem is desired which 
minimizes the quadratic function 

'Kp) = (y-Ap)"'' W (y-Ap) (A1.3) 


where W is an (m x m) weight matrix. In our applications, the weight matrix 
is diagonal with elements equal to the inverse square of the observation noise 
sigma, a , 


W = E[vvV^ = 
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A necessary and sufficient condition for a minimum of Equation (A1.3) is that 
its first variation be zero. 


6J(p) = 0 

This leads directly to the system of normal equations 

a''’ WA p = a''’ Wy (A1.4) 

and to the least squares estimate 


p = (A^)"^ A^ Wy 


(A1.5) 


The matrix A^WA of dimension (n x n) is called the information matrix. If 
in addition to the linear system (Al.l) we have an a priori estimate of p 
and an a priori information matrix denoted by p^ and , respectively, 

then the normal equations become 


(aWaq) p = aV + Aq Pq 


so that 

si = i-i^ = (aWAq)-I A'^WCy-A^^] . (Al.fi) 

The set of measurements in our application consists of magnetic 
anomalies in the total field, AB , and anomalies in field components, , 
Fg , F^ at geographic positions i . We consider two different parameteri- 
zations of the anomaly field: (1) a set of ND dipoles of magnetization Mj 
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with their direction fixed along the main field, where 


Ml 

M2 




and (2) a set of ND dipoles with components (m , m. » ) where 

j j ^3 
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In the second case there are 3 x ND parameters. Assume for the moment that 
there are no regional constraints so that all dipoles are independent. The 
parameter state vector p then contains all dipole parameters. Denoting the 
calculated i— measurement as 




= a 


ik Pk 


(A1.7) 


(summation convention assumed), it is readily seen that the elements of the 
matrix A are 



3Fi 


3p. 


(A1.8) 


and represent the anomaly due to the j— source at the i — position for 
unit magnetization. The calculation of these source functions are described 
in Appendix A. 2. The elements of the information matrix 


m 9F„ 3F, 


T n £ £ 

A.j = (A WA + 


(A1.9) 


and 


A W(y-y^) 


m 3F 

y W — -{y-y ) 

C— i J 


(Al.lO) 


where 


^0 = ^Po 


(Al.ll) 


are calculated and accumulated after processing each measurement. 


17 



Business AMD Techmologic4L Systems, Imc 


Now consider that the set of dipoles specified by the state vector p 
are resident in a total of N regions Rj;j=i,N further let the vector 
Pj denote the parameterization of region Rj . If the option to force 
the dipoles along the main field direction is specified, Pj has a single 
element; otherwise, Pj has three independent components of magnetization. 

The total state vector now is P with element Pj and a total dimension of 
N (or 3*N) and should replace p in equations Al.l through Al.fi . The 
il!l measurement is now 


Fj = I il )P. (A1.12) 


where the summation of 
elements of the matrix 


)lj is over all dipoles within region Rj . The 
A in equations A1.8 through Al.lD now become 





(A1.13) 


Advantage is taken of the fact that A is a symmetric matrix so that 
only the upper triangular portion is accummul ated. A Cholesky decomposition 
method is used to obtain the inversion of the information matrix, . 


The estimate error covariance matrix is 


C = E[(^ - P)(£ - P)"^] = A-^ 


(A1.14) 


and the solution correlation matrix p is computed as 
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The uncertainty in the computed angle $ between the main field in 
region Rj and the magnetization vector Pj for region Rj may be 
obtained from the 3x3 submatrix Cj from C corresponding to region 
Rj . Let 




(A1.16) 


where 


B 


1 



cos $ 



) 


B 


2 



(A1.17) 


B 


3 




cos 'll 



) 


and F , Fp , F 0 , represent the main field magnitude and components 
in the center of region Rj . Then 


a! = B^ C. B 
J 


(A1.18) 


The uncertainty in the magnitude jPjj may be obtained from equation A1.18 
by replacing the vector B of equations A1.16 and A1.17 with 
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(A1.19) 
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APPENDIX A. 2 SOURCE FUNCTION DERIVATION 


In this section we derive the expressions for the anomaly components 
and the anomaly in the total field due to a dipole at the earth's surface. We 
use a spherical coordinate system (r, 9 , i)) , where r is radial distance 
out, e is colatitude, and <|) is longitude east. Let primed quantities 
refer to the location of a point dipole, unprimed quantities to an external 
position at which the magnetic field arising from the dipole is to be 
eval uated. 

The magnetic potential at (r , 0 , (*>) due to a dipole source at 
(r‘ , 9' , <t>') is 


V = -M' • V (1/t) 


(A2.1) 


where l is the distance betwen the source and the external point. 

M is the dipole moment, with components (m^ , m^ , m^) . l may be written 

+ r'^ - 2rr' cos 


where c is the central angle between the two positions. Then it is easy to 
show that 


V = {m^(rA - r') - m^ rB + m^ rC}/£^ 


(A2.2) 




where 


A = cos 9 cos 9' + sin 9 sin 9' cos - (|)' ) = cos c 
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B = cos 9 sin 0‘ - sin 9 cos 9' cos ((|) - <|)') 

C = sin 9 sin (()) - (ji' ) 

For future reference, write A1 = A , B1 = B , Cl = C . 

Then the anomaly field vector is 

F = _vv = - ^ ^ I V 

hr ’ r39 * rsin93(|)J 

We win need 93 ^} (A,B,C) , 

They are 

3 A 

Yq = -sin 9 cos 9' + cos 0 sin 9‘ cos ((j) - (|>') = A2 
= -sin 9 sin 9' - cos 9 cos 9‘ cos (<}) - (j)' ) = B2 
= cos 0 sin ((|) -(})') = C2 


3A 

sin93(|) 

3B 

sin63(j) 

3C 

sin93(j) 


= - sin 9' sin (<|> - <|>' ) = A3 


= cos 9' sin ((|) - (j)' ) = B3 


= cos ((j) -<()') = C3 


(A2.3) 


(A2.4) 


(A2.5) 
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Define the following quantities: 


ni = r - r' A1 
D2 = -r' 

D3 = -r' A3 


FI = rAl - r 
F2 = -rRl 
F3 = rCl 


F 


r 


3V 

3r 


= m^{3DlFl/t^-Al}/t^ + rTiQ{3DlF2/£^+Bl}/£^ + m^{3DlF3/£^-Cl}/il^ 


F 


3V 

e ■ r39 


= m^{3D2Fl/)l^-A2}/t^ + m g{ 3D2F2/t^+B2} + m^{ 3D2F3/jl^-C2}/t^ 


F = - 

(j) rsin93(j) 

= m^{3D3Fl/t^-A3}/)l^ + mg{3n3F2/t^+B3}/t^ + m^{3D3F3/t^-C3}/£^ 

We now have equations for the components in the form 

F, = + mg d^2 % ^13 

^9 ■ ’'’’r ^21 '’’9 ^22 "’9 ^23 

*^<{> ^ "'r *^31 "’9 ^32 '”<!) ^33 


(A2.6) 


(A2.7) 


(A2.8) 
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Then 


3F,. 

r _ . r . r . 

3tti^ " °11 arrig ~ ^^12 " °13 


! 1 ? d 
’ 21 


d 

Smg ■ 22 


d 

am^ 23 
9 


(A2.9) 


3F. 3F, 3F. 

i = d ^ = d ^ = d 

3m^ °31 3m- °32 3m^ ^^33 

“09 

The anomaly in the total field is 

AB = F sin I + F- cos I cos D + F.cos I sin D 
r y <|) 


where I and D are inclination and declination of the main field at the 
point of evaluations. 

Thus, 

3 AB 

— = d-- sin I + d_, cos I cos P + d-, cos I sin P 
3m„ 11 21 31 


3AB 
3m 9 


dj 2 sin I + d^j, cos I cos P + d ^2 cos I sin P 


(A2.1D) 


SAR 

= dj 2 sin I + d ^2 cos I cos D + d^^ cos I sin D 

The above partial derivatives are used to form the Jacobian matrix as 
described in Appendix A.l for the case in which inversion of field 
measurements to vector sources is being attempted. 
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The formulation is different for the case in which source directions 
are fixed a priori and source magnitudes only are solved for. In this case, 
we write 


= M(sin i) d^^ + M(cos i cos d) d^^ M(cos i sin d) d 


13 


Ff, = M(sin i) + M(cos i cos d) d^^ + M(cos i sin d) d 


22 


'23 (A2.11) 


F^ = M(sin i) d^^ + M(cos i cos d) d^^ + M(cos i sin d) d 


<t> 


32 


‘33 


where i and d are inclination and declination of the source. Then we form 


3F 

- d^^ sin i + d ^2 ^ ^13 '' ^ » (A2.12) 


and the other partials similarly. The anomaly in the total field is 


aB = F sin I + Fq cos .1 cos D + F, cos I sin D 
r 9 <p 


= M[sin l{(sin + (cos i cos d)d ^2 + (A2 13) 


+ cos I sin n{(sin i)d 2 j^ + (cos i cos d)dp^ + (cos i sin d)d^^} 


22 


23J 


+ cos I sin 0{(sin i)d 2 j + (cos i cos d)d ^2 + (cos i sin d)d^ 2 l^ 


32 


‘33 


from which we form 3AB/9M . 
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APPENDIX A. 3 PROGRAM INPUT 

Program input consists of a main field spherical harmonic model, an 
overall grid of dipole locations, data at the overall grid locations, 
parameters defining a selected sub-region of interest from the overall grid 
and the definition of the mosaic regions as subsets of the selected 
sub-region. The overall grid of dipole locations is a 6D x 32 array covering 
the U.S., and the selected sub-region of interest is identified by the 
parameters II, 12, J1 and J2. 


OVERALL GRID 



The program numbers the dipoles internally by preceding through the 
rectangular sub-region sequentially from left-to-right from bottom-to-top. 

This is the numbered grid by which the user must define the mosaic areas . The 
sub-region must be completely encompassed by mosaic regions, but a particular 
mosaic region need not be simply-connected. As an example, consider the 
16 X 16 sub-region shown in Figure A.3.1. The Region I mosaic is defined by a 
sequence of entries giving the beginning column number, the number of columns, 
and the row number for each contiguous row segment. For example. Region I is 
identified by 


26 




Bvsi\ess ai\d Technol(X}ical Systems, Imc 


Beginning No. of Row No. 

Column No. Columns 


1 

1 

1 

1 

1 

1 

1 

1 

2 

4 

9 

4 

9 


11 1 

10 2 

7 3 

4 4 

4 5 

4 6 

4 7 

7 8 

6 9 


3 

1 

3 

2 


10 

10 

11 

11 


The program will also expect the total number of entries for defining the 
mosaic region to be given. 
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Figure A. 3.1 Sample sub-region grid showing a mosaic 
region definition. 
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Program Control Variables 


For the most part, the software is controlled by variables defined 
internally in data statements. In the main program, the following variables 
are set: 


JOPT = 0 
1 


estimate source components 
estimate source magnitude only 


LI 



do not process ABp data 
process ABp data 


L2 



do not process AB 9 data 
process AB 9 data 


L3 = n 
1 


do not process AB<|, data 
process AB^ data 


L4 = n 
1 


do not process aB data 
process aB data 


Moreover, the variable NDIM is set which is the dimension of the upper 
symmetric portion of the normal matrix. Note that this value must be at least 
as large as required by the problem to be estimated and the array D must be 
dimensioned to at least this value. 


In Subroutine FUN, the following variables are set: 


II 


the grid number of the overall grid which defines the 
lower latitude of the selected sub-region. 


12 


the grid number of the overall grid which defines the 
upper latitude of the selected sub-region. 
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J1 the grid number of the overall grid which defines the 

western longitude of the selected sub-region. 

J2 the grid number of the overall grid which defines the 

eastern longitude of the selected sub-region. 

In Subroutine DATA, similar variables are set which define the 
sub-region over which data will be processed. Note that the data sub-region 
need not be the same as the dipole sub-region. 

Program Input Units 


Data input to the program is accomplished on units 5 and 9. The unit 
5 input is as follows: 

a) Main Field Model in FDG format 

b) Latitude and longitude locations of the overall grid (60 x 32) 

c) AB data on a (60 x 31) grid 

d) AX data on a (60 x 31) grid 

e) AY data on a (60 x 31) grid 

f) aZ data on a (60 x 31) grid 

The unit 9 input is the set of entries read in Subroutine BLKS defining the 
mosaic regions via a 313 format. The data consists of the following: 

a) The total number of Mosaic regions 

b) The total number of entries for Mosaic Region 1 

{the entries for Region 1 consisting of First Column Number, the 
number of columns, row numbers 

c) Same as b) for Region 2 


z) Same as above for the last mosaic region. 
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APPENDIX A.4 PROGRAM LISTING 
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OF POOR QUAUTY 


/*J33Pia:i I.I3S3S100 
// EZEC 0F0BTa,?4aHaZa£f 

//sDaaca.sisis oo • 

* XBLaua (1520) , C0?DP(1920) 


mSIIISS'd! »' -"“I*' »' 

nT^e.cTnu 0 ]S« 0 i,DMia 0 ); 0 P( 680 ),SDD( 3 ) 

Qx (c30) 


QiaEHSIOH 

cgaacH/c?/ ci(6a&). 

coaao B/3L3KS/ BPcttia, sDcata, laLaoa, aaaai. 


— — ••••W •Mf MW •• w r A M m A A m flU 0« a u fl r. 

COMttOH/DAZlST/IDAI.IOUJl, iCa «/BD« FTO* BPD 

coaaoa/aac/BiAipii ,ax6» (32) ,atlv (3jf .St - - 
coaaoa/BMc/p (630) .bhp (630) , 5pdpi i6a0) - 
• 0P0P4 (630) .2 (2) ,f a,fT;pi?,i£,ap, ' 


.5 ' 


^oioli '(ild^ i)f Dpi^(6^ 


Jtft lOOUJ (^) , 

CQaaOH/DBAX/O 

DATA LI. 1,2, L3. 14, JOPI.JCMSra/l, 1,1, 1,0,0/ 

EATA SiSDF,S ifisa ,SIGa|,SIGSP/ll o! bl oJa! oU.O/ 

Jt-iaxa)® (j-i)*aDxa -(j**2-j)/2 ♦ X 
aoia IS xaa auasxa of sLsasMTs xa the no&an, aAiaxx amo a 

H0la.HPABai1lpAlaMf/2“*“ of IHE A36AX 0 

IS IHE HOaBEB OF PABAHETEBS ESTIHATED 


aoia>300 

00 422 I«1,660 
CP (I) *0. 0 
422 COBTIBUE 

SEX EAIA COaPOBEBI BEXGHXS 

SHTDF* 1. 0/31 GOF 
■X0F3SaI0F«SiiX£F 
SaTBB>1.0/S2GBfi 
kXBBaSBI3fi«SaX££ 

SHXaXal. 0/SIGSI 
HIBT»SHI3X*SaIEX 
SaXBP« 1. 0/SIG3E 
BXBP«SUZ3P«SiiX£P 

SEX OP SOOBCE A&EAX 

CALL Fl»2 (JOPI,JCaSXB) 

I^SOLUTIOH IS FOB SOOBCE COaPOHBHXS, COHPOXE AHGLE 
SBTaZEH A PBIOEX ABO BAZB FXSLO 

IF (JOPZ.E0.1) GO TO 16 
00 15 I«1,BO 
II-I 

CALL ABGL(ABG,IX) 

COBXIBOE 
COBIIBOE 


15 

)6 


SEX OP OAXA SEX FOB INPOX 
CALL 0AXA2 


JCO BP B«0 
JCOBPB«1 


BO PBIBX FOB COBBELAXXON BAXBXX 
PBIHI COBBELAIIOH BAXaiX 


XgPX>0 XO ESXXaiXE SOOBCE COaPOBEBXS 
J0PX*1 XO ESXZaAXE SOOBCE aAGIZXOOES OHLI 


JE BSXB>0 
JEBSXB>1 


BO COBSZIAXBX 
COBSXBAIIEO 


L1>1 XO IBPOX BAXXIL FXSLO COaPOBEBX 

L2>1 XO XBPOX SOUXB FIELO COaPOBEBX 

U*1 XO XBPOX EAST FIELD CQBPUBEBX 

L4>1 XO XBPOX ABOaALX XB IHE XOXAL FIELD 

FIS BEXOBBS COBB ESPOBOXB6 COHPOXEO QOABXXXXES 

25 FOBBAX (4X5, XI C) 


icoaa,. 


,3P (32) 

0) , 


THE 
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OF POOR QUALITY 


L 

C 

c 

c 


COSIBACT DFDP /iCIOB A3 DiCTA lEO CCMSXBAiSXS 
OxOP IS ia£ ABBAZ Of PABXIAi^ 

IhE aOSAIC BBSXOUS OBFIMSD lii lEBBS OP XH£ IBPUX 3BI0 

i>X4p01*Ad 

94^ oX.KS(liDCttTa,SUa£I.,X£i.iiUilJ 
JlOCIIXa lo rafc ilUaB£& of SOSAIC aSGIOJIS 

MPCiixa IS ms of XMO£fBinj)A«aT PAaAa£X£Bs EsrifliTju 

MPCaXasaOCBTA 

IP(JOPT.ME.I) SPCBIE=3*lliX:NTii 

PfiXHT A PBIOEI COVASXASCS HAIBXX 
,CAIi COBLPS (D,Da,a?CSXS, 0 ) 

DO 3 J=l,JlDXa 
i D(J)>0.00 

NN^O 

BB^O 

XSBQ Bxsax BABE SIDE 
DU 2 J3l BP 
DU (J) >0. DO 
2 COilTXBOE 

DO 18 J>1,liP 
P (J) =0.D0 
CS (J)= 0 . 


13 

13 


CC^H^IHUE 
COHTXHOE 


CALE PQB OAIA PBOPIEE 

POSIIXOBS (LAIXTaBB,LO«GXTaOE,ELSVAIX01l)= (BEAT.iiEOil, 

6i*£Y} 


BHXrB(5,30) 

CAXi DATA 

DO 20 X>1,BPIS 

clxL (j&PT,JCliSZB) 

BB>BB > 1 

3 Foaa BESXDUALS (OBSEBPED-COaPaTED fXELD) 01 PUB PAHXOUS XBPOT CQBPS 


39 CAXi.DAiA '(IMIT, BPTS.iAST) 
I • dil 



DO 79 L*U4 

™ ” 

, DU 74 J-l.ap 
74 OPD?]J)«OPOP1 (J) «SBTBB 


7 2 COHIXBUE 


XP (X.2.EQ.0) 
NH«BHA1 
DO 7S J>1,EP 


GO TO 79 


75 0PDP1J)»DPDP2(J) *S«IBT 
DX-j; 3 TiX)-PT)*SHTaX 

GO TO 70 
73 COBTXHOE 

XF (L3.EQ.0j GO TO 79 
IIN«BHA1 

, DO 76 J»1,JIP 

76 DPOP ]j)>OPOP3 (J) «SHTBP 

77 COaiXliOE 

X? (L4.EQ.0) GO TO 79 

DO 73 J>1.BP 
_)*0FDPi. 


UU /« jai.BP 

78 DPDP]J)«D?DP4(J) ♦SHTDP 

/ n* /▼% ^ c>^iinf\ a 

70 COl 


COSIBAd OPBP VECTOR TO OIHSHSIOU BPCBTB AS DICTAXEB 31 COSSTBAXai 

al»i 
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OF POOR QUALil Y 


161 

162 


16J 

64 

165 


DO 165 j=i,acctia 
a2»auaijaEi. (Ji-i 
IPfJOPT .SB. \ 


COFDS 
DO 16 


>0. OD 

«ai ,52 




w « w I a— a i M a A 

CDFOa(J) *C2£r2 (J) 
CONTlaUE 


GO TO 162 


* OFDf (lELSUa (5) ) 


GO TO 164 
J3=3* (J-1) 
CDFDP JJ3*1 
CDFDF ‘ 
CDFOF 
DO 16 
a3»3* 

CDFDF 




tO.CDO 

•O.GDO 

•O.CDO 


.. a>ai,a 2 
(lELSua <ai-i) 

(J3*1)adj3tDP(J3»1J*OFDP (a3 + 1) 

CDFDF (J3«2} >CDiD£(J3>2WaF0P (M3^2) 


CDFDP (J3+3) 

cosiiilaB 
ai»a 2 *i 
ccHiiaaE 


•CDiDf (J3«’3)«OFDP (a3*’3) 


FDaa D» ( 1 H 2 afis tfcxor) asd o aAiaicES 

DO S J>1,SPCUTE 
D6 (JlaDaiJ) ♦CDFDP (J)*DX 

LC^LdC (J,J,NPCSTB) • 1 
DO 5 a>J,S£CaiB 
DC«DC ♦ 1 

5 ^ CDFDP (J) •CDFDP (B) 

79 couTiaaE 

, WXIE (6,7) Bfl.X (1) ,X (2) ,X(3) .SS(I) , 

IF IZAUX.Zu.O) GO TO 99 
IF (Ha.DE.gPdlZa) GO ZO 400 


3T (I),BP(I) ,DF(I) 


COMPOTE CHECK SOM COLOMII AMD IS7EET 0 aATBIX 

DO 6 (.xl.MPCHTS 
soao*u.D0 
DO 4 a*1,gPCSTB 
EC^LOC (L,a,liPCSTB} 

XF(L.L&.ar EOEOC (B,E,aPaiTB) 

soao> 3 aaD ♦ D(ec) 




SUMO 

TSIMT (MFC STB , HPCSTB ,D, DFDP) 


F3BM PABXaEIEB COSBECZIOM VECTOB DP FOB TdE ESTIMATED PABAaETEBS 

DO 530 J«1,MPCBTB 
OP (J) >0. 

DCdl 31 -» 0 . 

DO 536 Ss1,HPCliTB 

LC^LOC ( J , E , BPC raB)^ 

1 « EE. K) 


IF(J.EE.K) ECaLOC (K,J,BPCliTB) 
DP (J)aDP(J)+0(IC J*DB Sf 

DCC] J) -DCC (J) ♦ 0 EC) *OC (K) 


sa 


530 COMTIBOE 

COMPOTE ABO ?BIBT THE VECTOB C? OF ESTIMATED MOSAIC PABAMETEBS 
iBITE (6.58) 

FOBMAI (///• *•• SOEOTIOS ••••///• ' ,91, • PO* ,20X, 

**CP* -20X DP* ,18X, 'CHECK* , 18X, * TAB* ,20X, 'BATIO* //) 

DO 540 J-I.MPCBTB 


POaCP(J) 

CP(J)aCP (J) >DP(J) 

ECaEOC (J. J.MKHTB) 
7ABaOSuBt(E(EC)) 

^siTi^ivia/cp (J» 

UBIT£(6,53J[ J.M,CP(J) , 
FOBMAT (14,6020.8) 


53 

540 COMZIBOE 


DP( J) , OCC (J) , VAB, BATIO 


SIP AMO OB 0 B-COBTBACT MOSAIC PABAMETBB ( P ) VECTOB IBTO DIPOEE 
PABAMSZEB VECTOB THAT IMVEBSIOH IS OV^ SO BESIDOAES 


MAI BE CAECOEAZEO . 
ai»i 

DO 565 Ja1,B0CliTB 


34 



nnn r; 


Blsi\ess ai\d Technological Systems, Imi 


361 

5o2 


563 

564 

So5 

4 } 


.*I 2 =ilU»UaEl (J )-1 

IFIJOPT .Sc. 1 ) i 30 TO 562 

00 3oi a=ai,a 2 

ilELHOa Id) ) =P lIELSUd (d) ) +CP (J) 
CO NTIH OE 
GO TO 304 
DO --- 


0 565 a>d1,d2 
33=3* ti£iB Ja(H)-l) 
1 * 1 - 1 ) 


J3»3* 

P(H3+f) =P'(a3 + 1) *C? 

P(a3*2)»p(a3*2)+cp 
Pja 3 > 3 l»p{a 3 + 3 j+cp 
C0STIB03 
coaxiaoE 
ai»a 2 *i 
CO NIX HUE 

HaXT£(6.40) 

|OBflAT(/U,’llL OIPOiES KXTHXN THE aOSAIC HEGXOMS'/// 
•4X .'aXPOLE* ,3X,*PABH. ' ,3X, • fi EGIOd • ,3X- ‘LAT dX, 

LOH* , 5 X 4 «PiaA 5 £TEaV 6 i,'fio. '.sx.’fio.^sx.'So.*, 


* 28 X,*TAL 8 ESV/) 

PBXHI OUT VALUES JOE 
DO 


ALL DIPOLES ,iXXHXN THE BOSAXC EEGXOHS 


566 J>t,HO 

ULL aEGXOE{J,JOPT,Xa) 
IF (JOPl.Jv.il GO l5 5o7 
ia»(ia-n 73 -»i 
L» (J-1) *5 
LPLUS1^L>1 


567 


GO 

CONTXaUE 


LPLuii-tta^ *J»I-J?WS1.xa,ALAI(J),AL01i(J) ,P(L>1) 
J,LPLUS2,Xfi, ALAX(J) ,ALOH (J) ,P (L*2) 
J,LPLaS3«XB,ALAX(J),AL0li (J) ,P (L«3) 


«aXX£(6.60) 

LPLUS3«L«3 




6} 


?EBFOaa BOOEL SVALUAXXOa AGAINST OAXA 


CONTXHUE 

FOaaAT(X5,10F12.5) 

IHXXXALXZE VAEXaIxES TO P£ 

A-0. 

3-0. 

NA-0 
HB-0 
LAST-0 
HH-0 

iaxT»o 
DO 12 X-1,5 

i 4 i\i) - 0 . 

, SBEAN(Xl-0. 

12 SXG(Xl-0. 

COBPAEE OAXA HXXH STEXHBTXC FXSLO COaPUTEO FBOfl 


93 


PAEAaSlBB SOLUTXON 

BBXIE(6,90) 

Foaa^ (/iXf'coaPABE data axiu coapuiEO values fboh booelv/) 

IT* V / & B n 1 ' 


. iBXXE (6-80) 

ID^FOaaAI ('3* ,16X, •BADXAL',13X,'S0DTH* , 14X , • BASI* , 12X, TOTAL FIELD*/ 
ciLL.DAXA (XaXX«j(PTS.LASX) 


38 


DO 


300 1-1 
BLAX 
BLOB 
ELEV 



apxs 


CALL FOB iJOPX.JCBSZB) 

WjBB (X) ;fB,BI(D ,FT,BP(X) .FP.DF (X) , IC 


3«3*ABS 
BA-BA»1 
B3-BBO 
39 FOBHAI 


(OFlIl-XC) 

(SB (X) -Ffi) 4-ABS (SI (X) -it) «ABS |BP (X) >FP) 




I.2F7.2) 

(I. -rc)* ' 

T -FB)- 
,1 -FT)- 
BP U .-FPl- 




BB 

BT 

BP 


SB EAB (2) -SHE AB (2) * (BB (X) -FB 
SaXAB (3{ -sax AB (3) * (BT (X| -FTl 
SaEAB (4) -SBEAB (4) > BP (Xj -FP) 


XI -IC) -iTDF 
X -FB' 

I - FT 


-■XBB 
•aiBX 
-a IBP 


35 
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OF POOR QUALil Y 


GO 70 88 


83 0 COSTX^iUE 

17 iLA 8 X.£w» 0 ) 
A*A/?XOAT(iiA) 
8=B/FA0AT(NB) 

A«HTbf 
38 T 2 )s 8 A*iiT£B 
S»T 3)>tiA«*IB7 
SBI 4)s8A«iiTBP 
00 13 I»1,4 

SHEAS (S) »3HEAa (5i 
SHEAS (Ij ^SHEAH jll/SA 
S5T{5) »3HI (5) ♦iHT(I) 
IT (ailS 


SHEia'cS) 53 a£aii%f+laiAli\i) 


13 


SXa{l}<3uB. 
EflS(S)=sIaT 
SHEAS (S) ^SH 


saxiE 
HEXTE 
11 1 EOEHAT 
22 2 70 SHAT 
SBXTS ( 
70BHAX 
HSX7E 


33 3 


, — EAHJ3)/14*SA) 

SIG (5) *SCHI (BBS (5) ••2-SHEAH (5j ••2) 

S (D H 1 1 1 ] 1 

E ( 6 J 222 I E 


HEAS 0X7?Ea£SCE,SCAXAfi=',76.2/) 
BEAS 0X7?EB£SCE,VSCT0B-*,7 o.2/) 
(?aSlX).X=l,, 5 j 


44 4 
535 
66 $ 


‘^iSf‘t‘S 5 U 3 .'f‘ 5 f;.' 5 Ji 


fOEHAT SHEAS:*, 576. 2/ 

aHXXE (6,5551 JSXG (Il,X»1,5) 
70BHAX (‘O*,' £XG:^,3 ?o.2/) 
— -XE (5,66&1 (Q(I- '■ 

SHAIt'O* C: 

(AI.E^. 1 / sfiXE 


SBIXE 
70EHA 
17 


17 

X7 

X7 

X7 

X7 


:l 2 .eq. 1 
(E3.su. 1 

*14.E«. 1, _ 
J07T.EU. 1) 
.JOPT.Ew. 0) 






SBXTE (6,42) 
iSXXE (6,43) 
BIXTE i6,44) 
SBXXE (6,45) 
SBIXE (6,46) 


2 /) 


1 1 
10 
23 


SBXXE (6,11) SC 
7oaaAX ('3*,X5,* sooacss*/) 

UBIXE (6,10) si 
70BHAX ('0*,X5,* PABAHEXEBS* /) 

■BXXE (6,23) 

70BHAX ('0*,' PUSCa PASAHBTESS* /) 
iBXXE (6,95) (E (X) ,X>1,SP) 

aaxxE ( 7 J 32 J (E(x Ki=i'aP) 

35 70BHAX (378.2/) 

32 FOBHAX (7711.3) 

GOaPaXE PABAHEXEB SXASOABO DE7XAXX0S 
S>3 

X7 (JOPX.EU. 1) S*1 

00 22 J«1,S 

a«J 

soa-Q. 

00 63 X>a,SP,S 
S3 saa-saa^PiX) 

AVG>SUH/7E0AX(S0) 

saa«o. 

00 69 X-H.SP,S 
$9 SUa»SOa> (P (IJ-A7G) **2 

22 

■BXXE (6,64) {SDDiJ).,J«1,S) 

$4 FOBaAX ('3*,' PABAHETEB S0«* ,3E12.4/} 

X7 SOLOTIOS XS FOB SOOBCE COHPOSESXS, COapOXS ASGLE XS DEGBEES 
3EXHEES 7EGX3S SOOBCE OXBECXIOSS ASD HAXS FXELO OXBECXXOS 

IF^^JOP|.E^. 1) GO TO 67 

63 7OBaAT^{*0*,' ISGAE BSXBEES HAGS 7ECX0B ASO HAXS FXEEO*/) 
,,, «aiXS(6,11Q) 

113 70aaAIJ[/lX,'aiP0X£ lax. mm.* ,9A,*P1* ,6X, 'P2*,6X-'P3* ,6A, 
••P*,7X,^7ISC IISC*,5X,*FDBC tOEC* , 7X, • ASG PHI*,5x, 
•'SIGPHl* .SX.'SXGHAGV) 

00 65 I>1,sa 
IX-I 

CALX ASCL (ASG, XX) 

6 5 COSTISOE 
67 COSIISUE 

XF (JCOBPB.EQ.I)CAXX CQBLPB (0,OC,SPCSTB, 1 ) 

3 COSTISOE 


I 
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JSITE16.50) 

SiOIAi CCttfOSEKI'l ' 

J0f»J IM?UT 3O0IU COaPOiiEBT* 

• ciapaHEST' 

'J',* I3P0T A80a IH IP*) 

ItI'I soosce diGSitaoss oaLx*) 

•3*»* laVftfiT 10 souacE cuaPoaEais*) 


400 

50 

4 1 POSHAT 
'4 2 {OEEAI 
'4 3 FUfiaAI 
4 4 FOEiiAT 
43 FUBaAI 
45 poaaAT 
SIO? 

£30 

SOaBOaTIME IS£IT(LL,aa,A,B) 

= r5i"il'S“plS''if!§il§ic=iSIfiiJ‘' ■“*“ «0‘-K3p 

iDiaiSo^*'^ PHScisioa osi/,osaa,A2,i(i) ,a{1j 

i.TBOH>1 

17 (X.X..1.T. 1) GC TO 900 

Ia.1»U.-1 

K1*0 

La»aa-Ui 
iBO^oia 
DO 90 K«1,££ 

X30>I304>La 

ftPIT=HID*'l 

LE30«K-1 

X0£aiA(K?I7) 

DO 30 I«S,Ll. 
lao-xaoi'i 
osaa>o.oo 

2^130)30,30,10 

xxao«x-s 

00 20 X>1,LE30 


10 


20 CONTXSaS 
30 asua«A(XHD)-Dsoa 
I7(I.aE.KlG0 TC 70 
^ IF (OSUB) 90 0,900,(40 
40 coariHufe 

zoxa^AXOGio (lOL/saaL (osoa) ) -. 

X7(XOXa.£E.XOXGL)GO TO 50 

XDZGX«XOIG 

1TH0»»Z 

>0 D?xy«osoBT (osaai 

Al*n.JO/DPXV) 

A (xao) -> 0 Pxv 

B^K^»DfZ7 

38 ^tj®5P*“saa40Bi.2Ui)*osaa 

90 coaTzaOE 

00 152 K-1.XL 


iM I a~i,4 

DPX7*A(KPX71 
A1»( 1.00 /DPI., 

A2»i 1 . 00-0 Bi£ ( A 1 ) * DPI 7 ) /0PX7 


h). 


A(KPX7)»A2*0e£i i1) 
a(lI.-K*1)»A (KPX7) 
ziao-K-r 

XF(iBaO( 130,130,110 

no 00 120 l*i,££mc 

Xao-KPXT^Z 

120 ciOTzkE‘^^**‘““* ♦001,E(Al)*A(ZaD)) 

ia^Kf5y“’‘®° 

KPI7«KPX7-La-1-I 

LABF^XaO 

DO 151 Z»K,LX1 

uap.LAap-ia-x 

A (LAIp/* Al*l)Saa40BXiE (A 1) ‘DSD a 
Ht-7lSOp51,151,140 
140 DO 150 i«r,lEaD 
uao«zAiP4l 

150 ‘^«5g{,|*(WHO)>OSUa*AaHD>Z) 

151 coaixauE 

152 coarxauE 
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OF POOR 




OO lao 
LIi»D=KPIV- 1 
LAK?»iCPIV 
DO 170 l=i,LL 
nsm^o.DO 
DO 160 
LXSD>I.IliD«'1 

DS ail3DSUU*^A (L) *k (Lino) 

160 COtlTIHOE 

1(K?I7) »DS0a 
Lijio » Liiio«’La 
K?I7»KP1V>1 
170 CUUTmUE 
; siici -^ o.oo 

KPIV»KPI7*I.fl 
iHD-iJio^aa-K 
13 0 ccHTiaas 

iaii:2(6,921) IDIGL.LTBOH ^ 

921 POHSAT (//2X,* TSIMT ,2X, 217//) 

aZTOBH 

93 0 IDISL—I 
LTSOH>l 


raiT£i6.920) LIB02 
923 FOESAT (5 2, • • • • • 
STOP 1J 


IliyEBSIOti FAXI.20 AX BOW, 16) 


BEXOBtl 
END 

SUBBOilTXBE FOH (JOPT, JCHSTB) 

XOaPUTES TBB PABTIAL DESi»AXITES OF IflE OBSEBVAXIOMS 2BX AiJ. 
JIP(3I.ES IB Ai.A aOSAXC BEGXOBS AdO Xd£ FXEXD 7ALUES DUE XJ ALA 
DIPOLES BEPBESEEXED 3T THE rABAEBIEB 7ECXOHP 

SIDE SS (. " ■■ " " ' 

ESXOE PC 

B/DXPOLE 

coaaoE/cp/ CP 


DISEESIOE SS (226) ,CC (22o)^,CS (226) ,CXh(22o) ,SXH (226) 
DIEEESXOE PC01(3.3), A (3) 

COaaOE/OXPOLE/ ALAT (226 ) , aloe ( 22b) 


COaaOE/BLOZS/. S?CiiXB,E 0 CEIB,ISLEUE( 1920 ) ,EUaEL( 1920 ) 
COaaOE/SEB /0 ( 2 . 60 , 32 ) 

COaaOE^HC/P IsSo) . OFDP (6801 , QFDPI (6a0) , 0 P 0 P 2 ( 680 ) ,DFDP 3 (680) , 

* OFDP 4 ( 630 ) , 2 i ' rtFB , FT , FP > XC , EP , E 0 

• 3X3 (3,3) 
coaaoE/caix/DEOBaxM) 

COBaOE/EOIB/EDXa,ajl,EJ 2 
DATA APB/ 226 * 10 ./^,, ^ 

DATA APSXGB/ 226 « 1000 ./ 

DATA APSXCI/ 226 «. 373 / 

DATA A?SXGC/ 22 b*. 873 / 

LOC(i,J.aBia)-(J-i)*aDia -(J** 2 -J )/2 ♦ i 

3AXA**ixAili£SX SETS YALOBS FOB DIPOLE SUB-GEID EHXCH EECOaPASSBS 

ALL aOSAIC BEGXOES 

XI LOEEB LA2XXUDE GSXD EOBBEB 

12 OPPEB LAIXXUOE GEXO EOBBEB 

J1 HESXEBE BOSX LOBGIIOOE GBXO EOBBEB^ 

72 EASXEBE BOSX LOEUXXOOS GBXO EOBBEB 
3 1ST DXPOLE SPACXEG IE KXLOBEXEBS 
X3 0-1(2) 

XaaEE » X ( 3 ) , 

XHBE2 » TUBEE ♦ 6371.2 
CALL FLO (OEE,XBO,XBBBE) 

DC-OBE 

BC-XBO 

BC-XHBEE 

T-B>I(3) 

XCa«{90.-2 (1) ) *ABC 
CXC-COS (XCi) 

SXC-SXE(XCB) 

PB-0. 

FX-0. 

FP- 0 . 

DO*t> X»1,BD 
L- (X-l)*3 
CXDiSTH(X) 

SXO-CXfl]X) 

OPE- (X(2)*ALOE (X) ) «ABC 
C0P-<:0S(bPB) 

SOP-SXH(OPH) 
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OE POOR QUau j ^ 


OD-SS (I) 

3D-CC (I 
CD»CS]li 

Al=CTO*CTC*3ia*STC*C03 
5 1 aC IC*S TD-S IC *CT3 *CD P 
C1»STC*S0P 

*2»-SIC*CrD>CT.C*S2J*C0P 

B2»-SIC*STi)-CIC*CTD*CBP 

C2»CTC*S0? 

AJ»-STD*SDP 

B3=CTD«SDP 

CJ»CDP 

B»H*E-»T*T-2.*T*a*Al 

imr 

C=VOI./C 

P1aT*41-a 

F2— T*31 
PB»T*C1 

oi»r-8*ai 

02«-B*A2 

D3»-a*A3 

D1-3.*D1 

D2«3.*D2 

D3«3.*DJ 

B 1 - F 1 /E 

P2«?2/E 

F3»P3/S 

01 1- C*(01*F1-A1) 

012- C»fDl*F2^aT 
013»C*{01*F3-C1 
02 1*C«{D2*F1-A2 
022»C«{D2*F2*B2 
02 |»C*{o2*f 3-C2 
031»C*(D3*F1-A3 
032«C* t03*F2+B3 
O33-C*(03*F3-C3 

^ IP (JOPT.su. 1) GO TO 30 

J FOR SSTIBATliiJ SOOBCS CORPONEBTS 


0F0P1 (L»1 
DP0P1 L>2 
0 FDP 1 i > 3 , 
OFOP2 < L* 


■ Oil 
■012 

■ 013 
■021 


30 


0P0P2 1*2 »022 
OFOP2 L*3 >023 
0FDP3 X*1 =0i1 
07DP3 L»2 >d32 
DF0P3 L+3 >033 

0P0P4 1 >01 1*OC*021*BC+031 *EC 

0FDP4 L*2 >012*OC*022*SC*032*ic 
DF0P4 (L»3{ >013*OC*O23*aC+O33>EC 
SS*CS*S3 JI® ♦012*P IL*2) ♦D13*P (i*3) 

FT>FT*'021*? (L*1J*022*P (X»2i *023*P |l*3) 
FP»FP*D31*p]l»1)*D32*pJl*2J >033*1 1l*3) 

GO*fo°31 >07034 iX>2) *P |x>2j >0F0P4 (i>3) *P(L>3) >IC 

COBTIBIIE 


Z FOB ESTXBATZS3 SOUICS RAGBITOOBS 

>00*01 1+aO*D12>C0*01 3 
>00*02 1>30*022*CO*023 
I >00*03 1>50*032>CO*03 3 


31 


0F0P1 , 
0FOP2 I 
0FDP3 ( 
0F0P4 
FH>Pi 
FT>P I 
FP>? I 

rc«p i 

COBT3 


4 >00*03 1>BO*032>CO*03 3 

■ioFDPUlf Jja* >flC*0F0P2(I) >BC*0F0P3 (I) 


l*OFOP2 ;1 >FT 
)*DFDP3 i *FP 
|*0FDP4 1 >YC 


I ‘i^aiAX.SOUTH.AliO EAST 

c 12 IS ABCBAXI IB T£E TOTAL FIELD 


37 COBTIBUE 
BE TUBE 

EETBZ F0B2 ( JOPT, JCISZB) 
BJ1>J1 
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t '**“s 

OF POOR QUALaTY 


aj^-J 2 

CiiEaO. 

XMOO. 

xaa££=:o. 

luai^ii - IHa£Z * 6371.2 
CAI^ PLD (CS£,IHO,XHfi££) 

FLO TAK2S (OHEflSiO ,IH2£E) a (i.ATlTUDi;,LOKaiXaOa, ^^VAIImNI 

aaxjsjs loat.XiiO.iEHSET* (sia i.cos i*cos o.cgs i.sia oi co«foi20 

FBOfl A FIELO BOOEl, UHSiui (I,D) ABE (laCLlaATZOa .OECAiaAIIOB) 

BEAD SOURCE LOCAXICUS FOB FULL 60X32) EEGIOH laXO ABBA7 0(2, oO, 32 

00 1 11 J = 1 ,JT0T 

00 51 I = 1,6 

1 STB! » 1 ♦ 6*(I-1) 

XEMD * I3IEI ♦ 5 

Xiao - axao (xiao.xxox) 

51 coaxiafil* '“t 2 ,J,x),Kaxsxax,x 260 ) 

111 coaxxaul 

3 FUBflAX(3F10.S,2X5) 

2 F0BBAX(12F6.2) 

XSEAX SUBSET OF SOOBCES HBXCd BXLL CCnfBXSE XOXAL OF SOSAXC BEGXOas 
L*0 

OO 3 X>11,X2 
DO 3 J«J1,J2 
LaL*1 

OHE*0(1,J,I) 

T30>0(2« J»X) 

THBEE»0. 

ux,Loa souBCE eosiTXoas 


IF(L.EQ. Ill oa£>oaE * 
ir L.3Q. Ill XUO>TBO *■ 
12 i.Bw. 112 “ ‘ 

XFiL.EW. 112; 

12 <..Ew» 12 < 

ALAT (X) «0a2 
ALOa jL^«XHO 


oaE^oas - 

XHO>XiO * 
oas«oa£ - 
ZHO«Tao * 


1 

,45 
, 4 
, 4 
, 4 
> 5 


xa>oi 


»liC 


CXB(L) >C0S (XB) 

SXH(L) »sxa (TH) 

10 211 SOUBCE 7ECIQB OBIEaXAXXOB la dAXa FIELD OXBECXXOa 

XBBES » XiiBEE * 6371.2 
CALL FLO (QBE, IBO, XBBES) 

SS (L)>OaE 
CC (L) >XUO 
CS (L) >XBBE£ 

a coBTiauE 
ao«L 
ap»ao 

XF (JOPT.as. 1) aB>ao*3 
TOL>01SX*0XSX«40. 

xsxzs>af«(ip * i)/2 
xfIxsize.gz. aoia! stop 15 

ZEBC aOBBAL aAISZZ ABO A FBXOBX FABABETEB VALUES 

00 5 j«i,ap 

> ccaxiauE 

00 7 j-1 ,acxa 
7 oaoBax(j) «o.oo 

XFUcasii .Ea.a)BsiUBa 
IFiJOPI. Ew.O) so XC 100 

IF SXAXISXXCAL A PfilOBX COMSTBAXalS ABE XO BE XBPQSEO, PROCEED 

SEX A PBXOPZ PABABETBB VALUES ABD BOBBAL dAXBZX 

SI-1 

00 SO J»1.a0CHXB 
a 2 *ai>HUB&L (j)-i 
c? (j)>APa(j) 

LC>LOC (J ,J,BQCKXfi) 

oaoBax{LC)«oao£ax(LC) * i.ooo/apszso (J) ««2 
00 no a»ai,a 2 
p(i 2 LaoB(a])>APa(j) 
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4 } 

33 

’103 


113 


123 


C 

'•m 

c 


COSIISUE 

ai=d2+1 

CONTIiiaS 


BETUttM 

SEX A PiiXOal PAEAiSEXEa 7AEUES AA'i) HAIBIX 

L “-2 

ai»i 

00 153 X^I.HOCNTfi 
a2=ai*iiaaEL (J)-i 
Fi>a.oo 

F 2 = 0.00 

F3*3.03 

DO 110 a*Hi.a2 
F1»F1»SS (ISLSUaU) 

F2»F2>CC(XELaUBlaj , 

F 3 = F 34 CS Usi'NOilUl ! 

COMTIHOE 
Fi»f VauaEX ( J ) 

F 2 =P 2 /auaEL I J) 

F3>P3/auaEX (J) 


J3»3* (J-1) 

oo 12 U a>ai,a2 

a 3 a 3 * iXEiaaa ( ai - i ) 
p { a 3 * i)»APa ( J ) »Fi 
p(a342)«APa(j) *F2 


Pia 30 l.»APH 

cpariaoE 


J) *F3 




_ , ■ A 2 fl(J 

CP (J3*’2} >4PB (J 
CP (J343) -APB (J 
1»J3*1 

HF2*F2*F2 


• FI 
•F2 

• F3 

♦ F3*F3 


aF»0SQiI jfiI2]i 
- 1, ir»FT/APSXGB 


F2/APS1GB<1) 
F3/APSXGB jxi 


3I2J3 •- Fl * F 3 / aF / Oa / APSIGlU 5 
III 2 •- ?|/ H £ 2 > 2 ?/ fS|IGD g j 



903 

153 



OaOBBXlXO 1] >013 (1,2) 
oaoBax ixc»2 ) •exb (i .3) 

XP1»i ♦ 1 

XC>L0C (Li 1 , LPT, aPCHXil) 

OaOBaA(LC)>3IB (2,2)^ 

OBOBBA(LC * 1) •3TB (2,3) 

LP2*L ♦ 2 

LC>LOC (LP2,LP2,liPCnB) 

OBOBBX (LC) «STB (3,3) 

B1«B2*'1 
BBXX2 
BBXIE 
UfiXXE 
FOBBA. 

coax 

BEXOBB 

'OBPOIE AHGLE 6ETBE2B OXPOLE AND aAXB EIELO AMD THE STAaOABD OEVXAIIU 
EHTBX AMGL (AMfi,XX) 

LL^O 

X(2(>AL0MU) 

CAiL " FM ^[ lio !- ofx ] l )^ X ( 2 ), X ( 3 ), 1968 ., 50 , LL , A 1 , A 2 , A 3 , A 4 ) 

PUXC>AIAM2 (A2^-Af) 

C2--A1/A4 

C3«A2/A4 

C1>-A3/A4 

XDEC»PHI C/0. 01 74533 
TXNC«AfiSXM (Cl) /O. 0174533 
AMG*0. 

XF (J0PX.2Q. 1)GC TO 9 
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si 


OF POOR QUALjTY 


pi*p (i+i) 

B 1 *P 1 /EM 
B 2 »P 2 /aa 
BB=P3/Ba 


23 


14 

12 


1 j 


AMa»AacOS (31ACUS2*C2*B3*C3> /O. 0174533 
flSC^ABSIiJ (fli) /O. 0 174533 
|0|C»ATAS2 {S3-£2)/0. 01 74533 

"’‘lALPSEGioitl^JOPT^iif 

IP3»IP>2 ' 

DO 20 IC»I£,IP3 
DO 20 JC=IP,XJ3 
LCaLOCaC. JO, aPCHTfi) 
imc.£|.jc} ^ic=Loc(Jc,ic,apcjiTS) 

KO? (IC-iP+1,JC-IP»1)=DaofiBX(LC) 
COS?»31 •C1+52*C2+33*C3 
SIHP»Sw3I (1.0-COSP*COSP) 

A <1)*(wl-COSP*B1J /tfia*SXSP) 

A (2j » IC2-C03P*B2) / [£a*SISP) 

A ]3i ajC3-C0SP*B3j / (Ba^SXHP) 

SXSPHX>0. DO 
DO 12 L«1,3 
DO 14 K»l'3 

^Oj|^GP|X=SIG£HX>A(K) •PCOy(A,I.).*AiI.) 

SIG?HX*SaaT (SIGPHI)/ 1.01745 33) 

All) aB 1 
A (2) >B2 
AJ 3 Ub 3 
SIGIIAa>0. 

DO tS L>1,3 
DO 16 K»ll3 

SIGaAG^SISHAG * A (K) *PCOV(X,X) «A (X) 
SIGHAG-SiiSl (SIG8AG) 

GO TO 10 
com HOE 
BH-PjII) 

PI— ia*ci 
P2— Ea*C2 
P3» aa*c3 
PISC»XXHC 
POEC«TDEC 
C0H7XHUB 


10 


**PDK Toii »ALOH(I),Pl,a2.P3,aa,PIHC,IXHC, 


(?1,P2,P3. 

:ALAi,iIoa 

TOEC,IIHC 


■aAGBEIXZAXIOH ?SCtOa COaPOHEHTS 
>aA6HETXZATX0H 7ECTOS PQ3XTX0H3 
>D|CUllAIXOH,XHC^HATIOM OP EAXH FIELD 


1 

3 

35 


OF aAGHETXZAIXOM VECXOB 
COdPOMSHTS IF SOXUZXOH 15 FOB VECXOB 
MBGAIXfE aAGHSXXZAIXOH XF SOLUXXOH 

X;> FOB aAGHIXODE OKXX 

HG X5 AHGLE B£XB|£H HAXH FIELD (POGO02721 AID SQLUTXOH VECTOBS 
(19!*i*^F6.1 ,3X,4Fa.i,3i,2F6.i,3X,2F6u0,3X,F10.0, 


cuaaAi; i' u ■ 

.F 10 . 2 , 2 x‘f 10 : 2 /f 

?oaaix(XS.2F6.i, 

BETOBH 
EHO 


4F7.. 2 ,2f 6. a ,2F6. 0,F 10..0) 


SUBBOaXXHE iaOXX(B,aX3) 
CAxcaxAXEs xfiE aAzaxx bib* b x 


ST 


DO 10 J«1(.3 


13 

23 


“ciaflJol^®*^**^’ * B(K.X) •B{E,J) 

CUHTXHUE 

BEXOBH 

EHD 

SOBBOaXXEE FLO (OHE.XiO.IBBES) 

DATA L/1/ 

CAU FOG (1,0,0, CHE.XaO.IBBEE. 196B., 50,L, A 1, A2, A3 , A4) 

CHE"-A3/A4 

I¥0»-A1/A4 
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TaBSS«-»i 2 /A 4 
i »0 
SETaBui 
BUO 

SUdBOUTISB CATA {ISXI,APTS,;.A 3 I} 

SELECT DATA XB PBGFILES 07 ££ THB SDG-GBIi 3 SEdlUill OBBIBBO BY 


X 1 ,I 2 ,J 1 ,J 2 IB BAZA SZAZBdBHZ 
oXiiEBBIUM BE (3 . 60 , 32 ) ,01 ( 3 ,b 0 , 32 ) ,02 l 3 ,o 0 , 32 ) ,03 ( 3 . 60 , 32 ) 
UXBEBSX 03 L£l 3 ).aii 3 ) 

coflflos/saa/o ( 2 .So,i 2 [ 

COiUlOS/BEC/BXAi ( 32 ) . kx .06 ( 32 ) ,Ei.£V ( 32 ) , 0 f ( 32 ) , BB ( 32 ) .61 ( 32 ) , BP ( 32 ) 
COBHOM/CATASI/iaAZ.IwOa.fc.FA.FT.i'P 


BEAL «3 SZO.OHEAN.FI 
DATA XIOT, JIQT , X 1 , 12 , J 1 .J 2 / 3 2 , M, 06 , 22 , 32 , 43 / 


XT (lail.Hi 
XMIT -1 


. 0 ) GO TO 30 


£Q. 0 | 

aw.O) 


LL >0 
JJ »0 
K -1 
XT 
IT 

3*0 

50 COHTZBUB 
0*0 

00 16 1 * 11 , X 2 
I,*L *1 
U.*LL» 1 


' (iLB(K). 
,|L£ (k| . 


K»A ^1 

K»K *1 


Bi.AT(l)*D( 1 ,J,X) 
BLCB(XU 0 ( 2 , 3 , 2 ) 
EX.ET (L)*BZ(K) 

8 B a)*- 03 (K,J,X) 
3 T(t) •-OlfK.J.X) 

apiLi-ox (s.J.l) 


IS OFlLj-DElS.J.xJ 
3 U *3 J ♦ 1 
HPTS*L 

BBXTE^jS.a), OJ.BPTS.K 


PBOTXX.B' , 15 ,', 
' GO TO 30 


3*«1 ♦ 1 

poaaAi ('a*.* i 
XT (U..Ea>BxOT) 

XT (J.LE.J2) GO TO 17 
J»JT 

XT (LB (K).EQ.U) K»K *1 
17 COUTimJE 
aETOBtl 

20 CONTXBOE 
XABT*1 
aaiTB ( 6 , 21 ) 

21 TOaMAI i'O*, • EBD OF DATA SSIV) 

Li.iiia 


M 5 , ' POIBTS, XIEB*U 2 /) 


asiiE ( 6 , 22 ) azoz.JJ 

22 TOBaAI (*a*,X5,' POIBIS TB0a*,X5,» aaOFXLES* /) 
BEXOBH 

SalBZ 0 AZA 2 


liOIS«Q HO HOXSa 
ilOXS*1 HOXSS 
HOXS*0 
H 0 X 3 * 1 
aaBB «3 
K «2 

LE ( 1)«0 
LB ( 2) « 1 
LE ( 3^0 

DATA POXBT LOCATXOI GBXO 


03 a DATA POX ax SUE«GBXO FOB XHPOT 

SiaOLAXBO DATA BAY 65 XBPUT FROa THREE POSSIBLE ELETAZXOHS, 6 T(X)* 3 S 0 , 
450 , ABO SSO IB, OEPEBDXBG OB BHEZBEB L 6 (X )*1 
BIOT«BaB& OF LET ELS OSEO 
xrOX-IOIAL BHBB XBSOX DATA POIBXS 

XT»I 2 -X 1>1 

JX»J 2 -J 1+1 

aXll)» 35 Q. 

HI ( 2 ) • 325 - 

ax ( 3 ) * 550 . 

BTOT -0 
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DO 7 1=1 ,J 

IF JLB(I) .Ei. 1) BIOI=aiOT«-l 
7 COSIiaOE 

aior»ii*jr«sioT 

aSlD IH ALL 'DATA* TED) 

OB lE^ABOKALr III ToiAL Fiito ' 

(D 1,02,03) AAE (KAOIAl.,iOilIH,£ASI) COafOaSBXS 

00 136 I » 1,31 

DO 13b J = 1,6 

ISTBT = 1 ♦ 11*(J-1) 
lEBD = ISIHT ♦ TO 
lEHO = aiaC {X£aD,30) 

|SA3^5j 130) (OS (K,JJ,I) , JJ*ISZBI,ZEaO) 

CCNIINOS 
DO 146 I * 1,31 

DO 145 J = 1,6 

ISTSI » 1 ♦ 1 1* (J-l) 
lESD = IStEI ♦ to 
iiao = aiao uEao,6Qi 

aHAO^S^UO) (El lK,JJ,i) ,JJ=XSTBI,IEMD) 

COMZiaOE 
DO 156 I » 1,31 
00 153 J = 1,6 

isiai « 1 ♦ il«fj-i) 
lEHD » israi ♦ to 
ixHO f axao iisao.60) 
ccaTilfis tD2 (i,ju,i) , jj=isiai,xE)io) 
ccaTxaul 

DO 166 X * 1,31 
DO 165 J « 1,6 
ISiaX * 1 ♦ 11* (J-l) 
laiiD = isxai ♦ 10 

XEBO « also (IEHO,60). 
coaiiaiE li.W.D , JJ“isiaT,XESD) 

coaTzaul 

poaaAxinF7.2) 

Ix (BOXS.EO.O) so ZO II 


135 

133 


145 

143 


155 

135 


165 

166 
133 


ADD 30X33 70 •DATA* 

OH3Aa«0. 

3X0*6. 

XI*231457 
DO 9 I«1,XI0I 
DO 9 J*1,JX0X 

CALi. GBOiaX (Xl,SID,Da3A3,FX) 
'• ,01 |k,5,X)>PI 




13 
1 1 


30 


- (xx.srb^ ^ 

3X0=1. 

U*80911 
DO 10 X»1,II0X 
DO 10 J*1,JX0X 

CAXX Guoaax <XX,3ID,DB£AN,FZ) 

aExuaa 

£30 

sJOBOoxiaE saoaax (i7,3Xo,oa£Aa,FZ) 
zaPLxcii a£Ai*e {A-a,o-2) 

a£AX.«4 I 
A=0.00 

DO 50 Z-1.12 
CALX BAIOO (XT,1X,I) 

X7*II 

A*A4I , 

PZ»(A-6. DO) *SIC*0B£Afl 

BBxaaa 

SHO 

3DBB0aXiaS CCBX?Bp,S,aOB,HX) 
a£Ai*3 0(1) ,3(1) .dopaia 
S=PXB III b I=BOS PXB J*COX PXB 

D AEaiZ HOLDS KOBaAL 2QUAXXOH3 OB COVABZANCB aAXBXX 
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- v” -■ 

OF POO^ QUfX:-'Y 


s IS psxsTuai laaAi 

Xf(HI.£w. C) GO 10 320 
DO 300 J>i,sca 
LC = iJ-1) •SOE- (J*J-3*J) /2-^ 
DO 300 ll*a,}iOi 
u:*LC*i 

330 D(LC)aD(i.C)/CSCHT(S(tt)*S(Jn 
■" CONTXH02 


323 


330 


K*1 

aail£(6,720) s 
KBIIE i[6,730) D (K) 

DO 400 J>2,l:0E 
S>0 

DO 3S0 X«1,J 
3(X)»oTS) 

X^iCt’iiOB-X 
•SITE (6,720) 0 
«aiTE(6,730) (£(X),X=1,J) 


A • 4 

43 0 COUTZaUS 
iETuaa 

720 roaHAX(/* J **,I5) 

730 FoaaAT(* •,13E10.2) 

EED 

saaaoaTXKE Bi.as(HOCBia,auaBi.,iELaaa) 
ozaBiiSioH baaExh920) . xEi.aufl(i92ai 
CO aaoa /B ox e/ BDxa , 0 1 , 0 2 

SUBSOUTIHZ Bi.ES IS USED TO ISPUI THE IHFOBBATIOu 
DEFXaXllG IH& OXPOXES COaEaiSIBG THE 90SAIC EEGIOBS 

saaaooTiNX ai.xs task is to oETaasxsA the 

EIEaSST SCBBEBS OF EACH COKSTBAXHT BUCK. 

COaBEBZlZ TU& iOBBOUTZliE IS GEAESO TO 
A GAXO OF OlFOIES 32 X OO IS SIZE. 

IllFDZ FOR THIS SUBBOOTIBB SHOULD BE AS FOLLOHS: 
lOH AHO COLOaK HUBfiSBS ARE SELATIVE 10 THE OXFOLE SUB>GaiO 
BEGIOO OSFZaBO IH SUSROUTXHS FUN 

BUBBEB OF HOSAXC BEGIOBS 

aOH BUBBEB OF FiaST BOSAXC BLOCK , BOBBER OF EBXBXES FOB BLOCK 
COL BOBBER OF FIRST EBTBI, BOBBER OF COLS IB FIRST SBXRI, 30a BOB 
COL BUdBER OF SECCHO SBTBI • BOBBER OF COLS IB SSCOBO EBTBZ, BOH 

• 

BOB BOBBER 3F SECOBO HOSAXC HLUCK , BOBBER OF EBiaXES FOR BLOCK 
COL BOBBER OF FIRST EHZBZ , BOBBER OF CO^Ci IB EBIRZ.ROB BOBBER 
COL BOBBER OF SECOBO EBTRZ , BUB OF COLS IB £BTaz,ROa BOBBER 


13 


BCL-J2-J1+1 
XTOTAL«0 
OO 10 1*1, 1920 
■OBELU) >0 
lELBUB (X)*0 
COBTXBOE 


BEAD TBE BUBBEB OF BOSAXC BEGXOBS 
BEAD (9,50) BDCBTR 


BRIZE (6. 65) BDCBTR 
6a FORBAT (//•*••*•***•••//' 
DO 40 I«1, BDCBTR 

READ (9,30) XfiOU.BOBROH 
XROH«XROB-1 


HO. OF RSGXOBS IS ',15///) 


OO 30 J«1,BOBROR 
IROB>IRQBAl 
READ(9.50) XCOL,BaBCOL,XBUa 
BUBEL (1) «B0B£L (I].«HUaCUL 
XSTHT»(lEOi-T7*BCL ♦ XCOL-1 
OO 20 K«1,iOBCOL 

lELBUB (XT0TAL>K) ^ISTRZi’K 
XCK»1C0L-1 *K 
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23 

33 

43 

63 

13 

43 


53 

55 


ITK=IT0XA1 ♦ K 
ISK*IS2ST ♦ K 
COiiXIJIUE 

IX0rAL=IX0IAi4NU«C0I. 

couxiHas 

CO MX! HUE 

EOEUAI (102., 6113) 

HHIXE (6.15) 

FoaaAi (// 132 , •EEdxoii • , 132 , • [jaaasa o? di^oees’//) 

DO 45 IT^I.UCCfiXE 
ilhlXE (6,55) II ,»ua£li(11) 

COSTXaUE 

REXUBE 
FOBaAI( 3 I 3 ) 

FOBiliX (52. X9, 62,115) 

luaaouiXEE fdg (j.aa.aEXi.DLAi, DtOHG.g.ia, aM 2 ,i., 2 , 2 ,i,F) 

c*******^**** ***************************** ,****fT!l*-**l*.f 1 


c 

c 


J.EQ. 0 laPUXS lilXIUDE 6 JaALIIXaDE (Ka) RELAIX7E XO aELXPSOXD 
(GSODEXIC COOROIEAXES) 

J.SU. 3 03I20T FIELD COaPOSEETS SORIH.EASI, VESTXCAL 

IB GEOOIIXC COOaSXEAXSS 

J.ME.O LAX.SLOliG XE SPUERXCAL COORDISATES, 3 = GE0CEBT£XC RADIUS (KJ) 
1.31.0 DUIPOX FLO coapoatais Muaia.EASI, VtRIXCAL xa aPaEBXCAL COOR 

la.SU.O OSE DEFAULT 7ALUES AE*6373.16,FLAI»29a.25 

la.a&.o xapai talues for ae.flai or £xB.st call io fdg. 

«E2X. Ew.o DO BOX READ laPUT VALUES FOR EXXERBAL FIELD PARAREXERS 

rheh l is greater XHAB 0 

iEXX. EQ.O DO BOX EVALUATE E 2 TERSAL FIELD FRCa 80DEL „„„„„ 

BE2X.RE.3 READ IKPUX VALUES FOR EXXERBAL FIELD PARAREXERS aa&a 
L GREATER 0 

a EXT. EE. 0 EVALUATE EXXEREAL FIELD aOOEL 


C 

c 

c 

c 


C DLAX 


C 

c 

c 

c 


c 

c 

c 

c 


3LOEG 

3 

iBA2 

iEAXX 

saAixx 

Baxxxx 

!( . EQ. 0 
(.BE .0 

rzsBo 

IRAR 


(KE) a HER 3^0 
UE) RHEE JsI 


GEDOEIIC LATITUDE IB DEGREES aHEE J^O 
GEGCSHXRIC LATITUDE IE DEGREES UREE 3=1 

LOBGITUOE la DEGREES 

GEODETIC ALTITUDE 
GEDCEEIRIC RADIUS 

BAXIHUE DEGREE ABO OHOBH OF COBSTAET XERBS OP FIELD aODBL 
a " FIRST ORDER TIRE " ! I 

N ■ a SKCOBD a a a a a 

a a a lUIRO a a ■ a a 

FIELD BODEL COBFFICIEHTS SCEBIDT BORBALIZED 
FIELD aCDEL COEFFICIBBTS GAUSS aORBALIZED 

EPOCH TIBS FOB FIELD BODEL CQEFFICIEIXS 

BEAR RADIUS USED IB FIELD BODEL POTBBTXAL EXPABSIUB 
(DEFAULT = 6371.2) 


lODEZT.EQ.O BO EZTSRBAL FIELD SOLVED HITH BODEL 
SSIO EXTERBAL FIELD SOLVED iillH BODEL 


lODEXX 
L.EO.O 
L.GT 


6 . LE< 


lEAD^li^FIELD^aODBL AHD EVALUATE FIELD 
EVALD AXF FIELD AX OLD TIBE 


t 0^0000 000000 00 00 0000000000 0000000 ^ 0 ^***^***^*******^******^^** 
EwOIViL£ilC2 (SBflXXp# 1UXG(1 ,1) ) 

^OlffifisifiB^fiHsJfai fox (18,18) ^HfllT (13, 18) .AID (33) 

DIBEBSIOB GTTT (8.8) ,aXl(18,ll) 

DATA XFRSI/O/ 

CATA AE, flat/6378. 16,298.25/ 

DATA TLAST/O./. , 

DATA IAaAR/6371-V,,„ 

IFdFRST) 110,100,110 


EQUATORIAL EARTH RADIOS ABO FLATTEBiaG FACTOR 
USED IB GEOXEXIC'GBOCSRTBIC COORDIBATES. 

THE BODEL ITSELF IS IBOEPEliOEaT OF THOSE 


C 

c 

c 

c 

c 
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OF POOR 


100 IF (nfl.ms.O) a£££(5« 101) ae.fiat 
13 1 FoaaAT J.1 1^216, 1) 

139 FoiaA:^i^/5xl 'CCalTisTS USED : V.22I.'iiOUAXOaiAI. EABId EAOIUS 
3r8.3/,22X,'EAaia iSCIPHOCAL FLAIlAEIHS ',Fb.1//) 

IFEST»1 

FLAT^I. -1./II.AI 

£ 1 > 0 . 

E2-0. 

E3^0. 

U»AE**2 
A4aAE**4 
B2»iAE*PI.AT) **2 
A2Bi-A2* (l.-Jl.AT**2) 

A4 34* A4* i 1 . - f L AT**4) 


110 IP 


3 

103 

104 
103 


7 

106 


107 

102 

3 


Pil 17 19 17 

BE AO^ (i^i) lui , S haIt, am AATT XTTI , aOOEXT ,E , T2EH0 ,ABAB , 
Aia (I) -1*1.10) 
aa.t (4I2<2I2,2F6.1, 12A4.A2) 

(ASAa.E;.0.) A3Aa*tAsA£ 
iOiS.lOj) aADCI),I*14,33) 


'it 

IP 

bead (S. 103) 
FOEail (20A4) 
1*0 



k£f20l4/lo 



BSA.'SXZEBHAL FIELD 
SSX«'B0BBALI2 AXIOM (E* 

S5X, 'FIELD HO DEI EPOCH '.FE. 1/,^ , , 
SSx!'FIELD HOOBL HEAH BADIUS ',F6.1//) 
aAXM«0 

TEHP*0« 

BEAD J§«6) ■,fl.aiH,B)ia,SXMH,aZHH,GZXMH,aXT!lR 


FOBHAI (2I3.6P11.4) 
IP (H.Li.O) SOX07 
HEXM* (HAZO (M ,dAXB) ) 


ST ji,a)«STMB 
aiT(E,d) «3TTaa 
TEai«Iaizi (XEaz. AES iSXMa) } 
IF (H.EQ.1) SCXOS 


S(a-1 ,A) «HHB 
GT (H-1-il)*alBa 
GTT(H-1,M) *aiTl 


GO TO S 


axiia 

IP (iaXTTX.EQ..O) GO 10 107 
BEAD (5,6) M .a.SXXZEH.aTXXMa 

-i.O) GO 10 107 

‘ SXOP 106 
XXTIB 
GO TO 106 


IP(M.'EC.O) GO 10 107 
IPW.GT.8) 


(a,ai*GXXTIB 

_7(a.Sw. I) sc : 

GIXIia-1.M)»UXXXBH 


GI - - , . 

IF (a.Ew< 
Gxxi (a-1, 
so TO 106 
COaXIBOE 


IP mooEXI.be, 01 BEAD(S,102} B1,£2,E3 
F0BaAZ(6X,3P9.2] 


FOBBAlfl^ib H B.6ZlaG,iaxlHH,9X2H6T,9X2aHX.8X3aSXX, 
S8Z3HHXI,7X 406X11, 7X4HHXXX//) 


DO 12 a«2,BAXH 

DO 12 a-i,ii 
ai»a-i 

IP (B.EQ. 1) GOTO 10 
(M.GX. HHZ'' 


9 

10 


11 

12 

108 

13 


IP(ir.GX:HHiXIXl aBITE(6,9)H,a,Gia,H) ,Giai,B) , 

5GX lH,ai ,GI (HI, B) ,GTX(Sl,ai ,GXT(HI,M) 

IF VIlilTi}' 8&ixi (Af 21) 

SGX]ll,H)^ar (HI. S) ,GTX (a,H) ,sfX(HI, S) ,GXXX (fi,H( ,GXXX (HI, M) 
POfiaix(2l3;8Pl 1.4) 

SO XO T2 
COHXIMUE 

IPMTSX.MHXXXZ) HBITE(6,11)M,a.6(H,H),SIiM.a) , 

*IP XaZIXT ) BBIXE (6, 1 1) M ,a, G (M, H) ,GX (M ,H) , 

^FoSAflki”? A“\iX,P11.4,11X,F11.4,11X,P11.4) 

CO BXIHuE 

IP (HOOEXX- aE.01iBIIE(6,108) B1,E2,E3 

FOfiHAX (//5X, 3HEIXPL0, ,3F10.2) 

FOBBAX (1H1) 
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14 


13 


16 

17 


20 

C 


21 


IP (TEMP. 2u. 0.1 L=-1 
IP (K.iiE.O) 001017 
SaaiT (1, 1) =-1. 

00 IS 

SaaiT (S, 1| =SHflIT (6-1 , 1) »FLOAI(2*a-3) /flOitl (»-1) 

saaiT (i.E) >c. 

JJ*2 
OO IS 

saaiT (s,a) =saflir (s,a-i ) *sofiT (pooai ( ia-a+i) *jj)/ploat ( s*a-2> ) 
saaiT (a-i,M) »SHaiT (a, a) 
u* 1 

00 16 s>2,aAia 
00 16 a>i,a 

GJN.d) =3 (H,a) ♦£aaiT(M, 

GT ja.ai “Gi (M ,aj *saaiT 
GTT(fi,ai «GTT (6,ai •saali(«,a) 

IP (iiaxxix.a7.o.ASio.a.i.E.d[Gxrt (ii,a)- 
IP (a. £0.1) G0I016 

?.SHaiija-i, a) 

Gxr ' 

IF , 

coilxiaoE 
x*xa-xzaao 
OO 13 a-i,aAia 
OO 18 a>1,B 
XGX-O. 

XUX-0. 

IP (a. Ew.1) GO 10 270 

IP (H.Gi.aaxxxi) GO lo 210 


j3) 

X (s,a) 


Gxix (6. a) «saaix(6, a) 


Ha- V i i caii a) Ssaaif (a^ 1', u ) 
ci(a-i,Ii) *Gii ta-1,11) *saaix (a-1 ,ai 
■ (aaxxxx.GX.o.Aao. a.LE.d) aXii(a-i,a)=G‘ 

\ Uiw ▼ M ft ^ * 


ixxx (a-i,s) •SHdix (a-1 ,a) 


210 

220 

230 

270 

240 

250 

260 

18 

19 



XHX» (XHX ♦ Gii(a-i,an*x 
IP (H.GX.aaizx) GO xo 230 
iaz»(iGX ♦ G7(i,an«x 

xa x» {tax 4-6 X ( a- 1 , a) j • t 

XGX«TGZ«G(a,a) 

Tax*xax»Gia-i,x) 

XU (a,a)«xsz 
XG (a*i,B)>xax 
GO xo 18 
coaxiaoE 

IP (B.GX.aaXTTX) GO TO 240 
XGX*Gixi (a.a) *1 
IP (H.GX.aalxZZ) go xo 2S0 
XGX» (TGX+'GTX J»,ai ) *X 
IP (B.GX. aaAXX) GO XU 260 
XGX> (XGX^GX (B,B) )«X 
XGX> XGX^GItifa) 

XG (a.ai-TGi 

COllTISOE 

XLAST-Xa 

DIAXB«0LAT/S7. 295779SD0 
5IliLA>SXa(0LAZ£] 

BICa&>0L0HG/S7. 295779 SOU 

c?a«cos(aLoaG) 

spa>sia(aLOHG) 

IP (J.Ew.O) GO2O20 

W IS GEOCEBXaiC SAOIOS UHEa J>1 

a*Q 

ci^siaiA 
GO XO 21 
szai.A2*siaiA««z 

W IS GEOCEXIC AZtlXODE 9688 0*0 
XlX«Q 

CUS1.A2*1.-SIBI.A2 
02 V2* A2- A2B2 «SiaLA2 
02B-Sy8T (0822) 

'• ' ♦DE2) ♦82) ).**2 

I.A2>SI2LA2) 

A4B4«SIH1.A2) /DEa2) 



aa AX*aiii o ( Ha x, bazb) 
H2XXP*!iEXI 
CALI a AGP 
I»BP 
P»B 
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yi^Ln'V 


22 

23 


1 


2 

3 


4 


3 

a 

C 

C 

7 



161 

165 

173 


IF (J) 22,22,22 

X--3T 

2 »-dB 

BETUBN 

Td2li2?0i2S FIE13 TO uSOOiCTIC OIEaCTXON^ 
SIUD«SI6IE«SI-StiBT (COS1A2) «CZ 
C0Sa»Sv8I ( 1 . 0-£IHD**2) 

X»-3 T*CO SO -n B* SI Si) 

z»BX*siao-aa*ccso 

BETOBit 

EBO 

S0BB0UTIB6 BEfiF 


CQaBOIi /C3BFFS/a(1S,13l 
CO MHOH/P LOCO M/ ST , Cl , S P H, C 


^oiaEBSiOB P(i 8 ,iSi,bp tis.iaj .CoMSTcid, ia;,s?(i 8 ) ,cPn&) ,fbi 

(P(1^1» .Ew. 1.0) so TO 3 

DP (1,lT»0. 

SP C1I-0. 

CP |l)»1. 

00 2 B> 2 ,ia 

PH 

00 2 a« 1 .B 
?H {a)»a-i 

COBS? (B,a» -PIOIT ( (B-2) ••2- {B-1) **2) /FLOAT( (2«B-3) • (2*M-5) ) 




, . isil 
CP I2)«cpa 
30 4 a>3,BaAX 

SP (ai »sp i2) *cp ia-1 1 4CP (2) ^sp <a- 1) 

la-il-3pl2i*spU-ii 

AB«A0B**2 

BT«0. 

BP«0. 

BB*0. 

00 a B>2.iiaAZ 

AB«AOB*AB 

00 3 a*1 ,B 
IF (M-ai 6.5,6 


&■ — a| w#w#v 

P(H, S] »ST*P iH-1,il-1) 

DP (fi.ai»ST*0P(«-1,M-n 
SO to 7 


♦CT*Pta-1,J-1) 

p ( H, a) >cx«p ( a- 1 , a) -co nsz (h, a ) «p (b- 2, a> 

BOTE : CCBSZ(2,1)«D 

0P^B^a|^»cr*DP (B- i,a) -st*p (b- i,a)i-co«sx (a,a) *de ib-2,b) 


F la.EQ.i) so xo 9 
2BP»S (B,a) «C£JB)4S(a»1,B) *s 
p«ap- (b,b)*sp (b) ^ (a-i,B) 
0 xo 10 

X£BP«G(B,a)«CP<B) 
BX«SZ 4 X£aP« 0 £ (B,k)«AB 


IF 

IS 

BP 

SO 


•SP(B) 

«CP (a))«FB(B)«PAB 


BB«BB-ZEBP*FB (b/*PAB 
3P«BP/SX 
IF (BEET. SX.O) CALL BAXFLD 
B-^BZ (BX«3Z4BE*BP4BB«aB) 

BEXUBB 

EBO 

SOBBOUZIBE fiEGIOB(X,JOPT.XP) 

COaBOS/ 3 LaKS/BPCBXB,BOCBIB,XELBUB( 1920 ) ,aUBEL( 1920 ) 
ai»«i 

DO 165 J«1,B0CBZB 
a2«ai>aaaEL (J)>i 
JBEG-J 

DO 161 B«al,B 2 

IFCIELBOB CB) .EQ.I) SO ZO 170 
COaXIBOE 
B1-B241 
coax II oz 

STOP 199 
coazzEas 

IP-JBEG 

IFJJOPZ.EQ.O) XP«3* lJBEG-1) 4l 

BEZOBB 

EBO 

SOBBOaXXSE ZZTFLO 

COaaCB/FLOCOB/SX,CT,SPH,CPa, a,aaAX,BX,BP,Bfi, 3,ABAB,E1.E2,E3 

I1«B2*CPH4£3*SEH 

T2»E1*SZ-Z1*CX 

zi«Ei*cz«>ri*sx 
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3fl=Ba-T1 

aT=BT*T2 
BBTUBII 
EBC 

// S££C 01,IBK20a,a£GIOB. SO^IOaOK 
//G3. PTOoFOOl 00 03H-»f9#aG. tjacOEP (* 

// DO 0S8»IC0aH. IMlTEr" 

// DO DSa=JCDHa.IHVBiX.Arti 

//GO. FT09P001 00 DS«»rCCA3.CL_ 

30T2FIXS sixoJIjiSiLiiroixI^Dxspiaafi'lABlL* (, , .OUT) 
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APPENDIX B. FLAT EARTH PRISMATIC MODEL 

The flat earth software employs the prismatic model of Plouff 
( Geophysics , Vol . 41, pp. 727-739, 1976) and the linear least squares 
algorithm described in Appendix A.l to estimate the magnetization of the prism 
blocks. The program will operate in two modes: 

a) adjust source prism magnetization vector with the direction 
forced to lie in the main field direction. 

b) adjust both the magnitude and direction of the source prism 
magnetization vector. 

Data input to the software is any combination of ARp , AR 9 , aB^ or 
aB . The local coordinate system utilized is such that x is toward the 

A> A 

east, y is north and z is down. 

The anomaly data input is accomplished in Subroutine DATA via unit 14 
in the same manner identified in Appendix A. 3 for the Mosaic Dipole Program. 
The data sub-region is controlled by the variables II, 12, J1 and J2 
defined in Subroutine DATA. The order of input is aB , AX , AY , aZ where 
AX is the anomaly in the north direction, AY is the east direction and aZ 
downward in the conventional magnetics notation. 

The input defining the prism corner points is accomplished in 
Subroutine BLOCKS via unit 12. Information is provided in the following 
order: 

a) Number of prisms [format (I5)j 

b) Number of corners for prism I [format (IB)] 

{lat, long, flat earth X and Y position for each corner point 
[format (2F6.1, 2Fin.2)1 

c) Same as b) for prism II 
• 

z) Same as b) for last prism. 
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The following variables are set in the main program: 

ND the number of prisms to be used. 

ZF the height of the grid of anomaly data to be used in 

kilometers (note, negative upward). 

Z1,Z2 the depths of the prism bodies in kilometers. 

JOPT = 0 estimate source components Px > Py » Pz 

1 estimate source magnitude only, with direction along main 

field. 

LI = n do not process ABp data 

1 process ABp data 

L2 = 0 do not process AB 9 data 

1 process ABe data 

L3 = 0 do not process AR^ data 

1 process AB<j, data 

L4 = 0 do not process AB data 

1 process AB data 
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APPENDIX B.l SOURCE LISTING 
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OF 


PG 




//YrORHal JOu (F8002,377,2) , ?LA2. EAaXB, IiaS= 101 ,u0) ,N0TIPI=IC0ai»,i.i.A5i=i 
/*J7BPABa &aE3£-c£XCa 
// £X£C OfOAlH.PAJiasliii:; 

//3ISIN DO • 

DiaEHSIDM si (5) ,Sa£AS (5> ,SIG(5» ,aaS(5) ,SDD(3) 

DiasasiDN olisSi .omUi ,op(i4 .gcn4) .jci-iui 

COaaOM/OAr/Df j60,31) - ax(o0,3 1 ,flP I60,3l) ,|Bj50,31) 

COaaOH/B»3/E iT4) ,OFDP (14»4) , 17 , Y£,^.x ,FH , FI,FP,IC, SP/MD, Z1 ,22 
CO aacM/0 BAX/0 


coaaoa/saia/aoia 

C0aa0S/PG3/aPEXSP, XL*T (41^, ZLOH (4) 
a£AL*a 0,0?0P1,O7OP2,OP0P3,DcOP4,OfDP, 

DATA I1.l2,J1,a2/06,22,32,48/ 

DATA JcSaPi/ 1/,L1, L2,L3, £4, JOPT/1.1, 1,0, 1/ 
DATA 0£X., AS «,ZSa/. 524 34,14.7944 0.2. 3 1220/ 
LOCU.J, 40IBj» lJ-1)*aDia-(J**2-J)/2*I 

SO IS Ta£ HOBBEE 07 PHISBS 10 B£ ESTiaAXED 

ao>3 

ap>ao 

CZaXEB CSOaOXBATBS 07 PBISaS 


OB,OP,SUHO,OC,OCC 


HAT (1 

Yi.ua n 

ALAI (2 


= 42. 
*90. 
= 35. 


ILoa {2 =73 
HAT (3) =37. 

TliUa(3)=S5. , 

17 (JOPT. £*.0) BP»3*8P 

apaxsp=ap 

BDxa»(ap*(ap4.ii) /2 

00 1 j=i,aoia 

D(J)*0.00 

7ACtoa«7.5*25.4 

0HB=0. 

tao»o. 

TaaEE=o. ,,, , 

XH£BE»IHB££46371.2 
CALX 7L0 (OaS,T»0,THBSE) 

CALL BLOCKS 
CALL OAZA 
CALL ABGL 

jcoBPa«o ao paxax fob coBBELAixoa aAxaxx 
JcoapB=i paXBX coaa^Lixxoa aAiaxa 

JOPT=0 ro ESXXHAXE SOU3C3 COaPOaBBTS 
JOPX=1 73 ESXIBAXE SOUiCE aAGBIXUDES OMLY 

L1=1 X3 xapox EAOXAL fIBLO COaPOBEaX 

L2=1 X3 xapox SOUXB FIELD COaPOaEBT 

L3=1 X3 xapux EAST FIELD COaPOBEaX 

L4»1 X3 tapux AauaiLY xa xhe ioxal field 

FOB BEiaaas ccbbespoboihg coaPOXEO siOAaxxxzES 


A=0. 

3=0. 

KA=0 

SB=0 

LAST=0 

aa=o 

XHXI»0 
do 14 1=1,5 

ll^a *0. 

14 SIG^I)=0. 

RB=0 

ZERO aiGHI BABO SIDE 

DO 2 J=1,ap 
DB (J) =0. 03 
P{J)»0. 

2 CUaXIHOE 
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?r-i r. i'' v-7'' 


OF POOR QUALiTY 


CATA HiiLIGiiT 

Zf=-325. 

HZBT& 0? paisas 

Z1*0. 

Z2=«0. 

00 20 J-J1 , J2 
DO 20 i*n,I2 

•5£i*XSS-34.6n*FACTOi 
I? (DF (J-Il .£0.999.) GO TO 20 
IF»( (1-1) *DEi.>I3H-§.52)*P4CTOH 


CAU. FUN(JOPT) 

A»A*AaS(Df IJ,I1-XC) 

3=B+AflS(aa (J,IJ-Fi *433 (BI (J ,1) -FI) +A3S (3P (J,I)-FB), 


HA>lilA«1 

NBs3B>3 


mm- 


?aj *(3a (J ,x) -fa| 


sill :3 111: 

SSEAHIT) «Sd£AH (1) + (OF 
SaEAN (2( >3 BEAM 12) * (BE 0,1) -Fa' 
SaEAH (3) >Sd£<ta j3) + (3T (J,X) -FT 
SaEAB (4) >3a£A3 (4} + (BF J,X) -FP 


aa»aa*i 

FOEa BESZOOALS (OBSEBTEO-COaPUTEO FIELD) OX FOB VAfiXOOS XaPUT 
DO 79 L«1,4 

, GO TO (7i;72,73,77) ,1. 

71 COHTIHlJS 

IF (1.1.S(i.a) GO TO 79 

DI*(B« (J,Ii-?B) 

GO TO 70 
7 2 COKCiaOE 

IF lL2.Zu.li) GO TO 79 
aBaHH'fl 

73 CONIXaOE 

IF(L3.Ea.J) GO IC 79 
!ia»aN4i 

0I»(BP(J,I)-FP) 

GO TO 70 
77 CONIXaUE 

IF (L4.£a.O) GC TO 79 
aa>(iti>i 

0I^0F(J,X)-YC) 

70 COBTIHUE 

FOBa DM AiO D BAIBICES 
DO S JJ>1,liP 

OB (JJ)30B(JO)«OFOP(JJ,i.)«Dr 
LC^LOC (JO, JJ,B£) - 1 
DO 5 a«jj,sp 
LC»t.C ♦ 1 

5 * ofdp (JJ,i.) * 0FDP{a,i) 

79 COailHOE 

BBIIE (6,7) aa,aa,I, J,XF,ZF,FB,FT,FP,BB (J,X) ,BI (J,I) , 

* BPlJ.lI ,0F (J.XI 

,7 FOBaAI i4I5,4F10.2,F10.1,5Z,4F10.2/) 

20 COBTIHUE 

I? (BH.LT.BP) GO TO 400 
A»A/FI.UAT(BA) 

3>£/FL0AX(BB) 

DO IS I»1,4 

:a(5)-»u(x). 


aas(i) *saaf iG (li/HA) 

SaSAB (5) «sa£lll (5) 4SaEAH (I) 

1 5 111*1) =SG|T®lai5l(^^ (I) **2) 

-SaEAB (5) ««2) 


Basl . 
SaEAH 


BBIIE (5,333) (BaS(X) ,I>1,5) 
BBIIE (B.‘144n£BSAa(X) ,1-1,3) 
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BEITE (6,555) (SIG (I) , I»1 ,5) 
iiaiXE (6,5o6) JvU),X*1.6) 

iL2.£«. 1 iillZ (6,42) 
(L3.C2. 1 BfiXTE 
(L4.SJ. V iEITE 
(J 0 Pt. 24 > 1 ) 8BXT£ i6,45t 
WSXIE (6,46) 

• fills ( 3 , 11 ) ED 

SEXIS (6,10) HP 
iBXIS (6,23 


IF 

IF 

IF 

XF 

IF 

IF 


sails (6',95) (f (X) ,X*1,HP) 

COBPUTS CHECK SUH CQXnHH &HS IHTEBT D dAISXZ 

DO 6 1.31 (HP 
SOHOsD.OO 
DO 4 H»1 ,HP 
lAl^LOC (i.,H,HP). 

XF (L.XS.H) XC»LOC(B,L,HP) 

4 SOBO^SUaD * D(XC) 

3 oc(X)3Suao 

CXU TSXH7 (HP,HP,D,DFOP) 

FOaa PAHABETEE COEESCXXOH VECXOB OP 

DO 530 J«1,fiP 
DP (J)30. 

DCClJ)i=0. 

DO i30 K«1,HP 
LC* E0C(J,X.HP1 
IF(J.EE.E| XC-’ 

DP(J)3DP(J)* 


Loc 


530 


DCC J J) -OCC (J) 4C 
CO NX! HOE 


(XC) «US 

»ciu:)* 




IP) 


DC (K) 


HEXXE (6,5dl 

*®*Fp“^20X, • 5 p* 1 18X, *CHE&£*., 18 X, • 7AE* ,20X, ‘BAXIO* 
DO 540 J31,HP 


//) 


P03?(J) 

P(J)»P(J)*DP (J) 
XC^Loc (J.J.HP) 
7iE-DSvai{0 ■' 
DC (J) 3i5(LCr 
BAXXO«7AE/P 




33 

540 COIXXHOE 

A>a. 

3-0. 

HA-0 

IB-0 

XASX-0 

HH-0 

XHXI-0 

DO 12 X«1,5 

iliiH*(i|j»o. 


(XCi) 

f ^ r 


,OP(J) ,0CC(J) ,7AB,BAXX0 


12 SXG(X) 

; CDBPAEE DATA BIIH SIIIHEIIC FIELD CQflPUIED FBOB PAEABEIEa SOLOTIOH 
C 

3 0 ^BBAi^J'O^ (16X,*BAOIAL* ,131,* SOUTH' ,14X,' EAST* , 12X, 

♦•TOTAL tXELfl') 

DO 800 J-J1,J2 
00 800 I-I1,X2 

XF»{ (J-1) •DiL*’XSM-34.61) *PACTOa 


__»{(J-1)#L. 

XF (DF (J,ir.EQ.999.) GO TO 300 
* -.-if •l)EL4YSH-9l52)*FACI0a 


IF 

HI 


liiV 
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CAIi_FaiiJOPT^ ^ 

(DF(J 

(88 (J 

HB-eii43 


BBXXE(6,a9) 

A-A4 ABS(0F (J,X1*ZC) 
B-S4AES(BB {J,X}-FB) 

HA-HA* 1 


JjBE (J,X),FB,BT(J,X) ,FT,BP(J,X) ,FP,DF (J,X) ,IC 
*AflS (BT (J ,X) -FT) ♦ ABS (BP (J , X) -FP) 


FOBHAT 

u 


w(3. 

2(6) 


U.3: 




3,4 (5X,2F7.2)/), 
(2,I)-lc)*(DF(J,X) -IC 
J,X)-FE)*laH (J,xi -FH 
(El jIx -FI •(3T|j.x 5-FT 
Vi (4)* (BP (J'X)-FP)*{BP (J.I)-FP 
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QUALiTY 


BtS/yESS 4JVD TECmOLOGlCAL SYSTEMS, LSC 


SAEkH (1) ^SiiZkS (1 } * (OF (J,l) ' 


33 0 COUTIiluS 

i»4/Fi04r,_„ 
B>B/FLOklidia 
DO 13 X»l,4 


fh; 

p 1 


13 


Sis'd^isla^lvai/iiA) 
saBAii (5) (5) *sh 

"^AS UI »sa*AH lij /s* 
(ij^saar das (ij *• _ 

; (5) =Sv<aT u (5) / (4*a A) j 
SaZAA (Sr«Sa&Ai3 (5i/(4«liA) 

) »SQBI (Kai (5) ••2-sk. 


EAH(I) 

SaZAH jli/UA 

aar Lzai (ij *»2-sa2Aa (ij **2j 


SAB (S) **2l 


SIG(S)»5QBI I 
SailE ( 6 , 111 ) 

aEirij6,222) 

111 SOtEki {'O',' e£AB OIFFBEEBCE . SOALAa^ • ,F b .2/] 
222 FOoBAI I'O'-' EEAtii DIFFEaEBCB , V ACX08=' , Ft>.2/j 
anTn,^ (0,335) (BBS (I) ,1*1.5) 
i»3',' fiaSi',5F6,5/ 

(6,444) (SaZABd) ,1*1 ,5) 

(*3',» saSAMi • ,5 f6,5/5 
(fpd) J=i, • 


333 

44 4 


WRITE , 
FOSaAT 
WRITE . 
FoaaAT 
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03 6 


cKJoaAi r u* , 
IF (ll.BQ. 1. 

” *L 2 . 2 w. 1 


HUIE (of41) 

, , WillE 6,42) 

8RI3E 6,43) 
(I.4.Ea>1) HRIIE 6,44) 


(JOPt.EA* 1 ) 
(OOPT 0 ] 


■ BITS 
4 BITE 


BBITE (6,11) dC 
11 fORBAI dO',13,' soaacES'/) 
8BXI2 (6,10) as 

10 FoaaAt (‘O', 15,* 

aaiTE (6,23) 


23 FORBAT (^3',' FBXdT PASAaBTEBS' /) 

dBIIE (6,35) 

3 5 FORMAT '' 

32 FOBMAX 


(7F11-3) 


' PABAaZXEBS' /) 

_RT PASAMB' 

(I) ,1*1, HP) 


IF SOLUTZOH IS FOB SOURCE COaPOHSBTS, COaPOTS AMaiB IH 
OEGBEBS SEIWZEil VECTOR SOUBCB OXRECTXOHS AdO MAXH FIEX.0 
OXBSCTIOH 


1) GO TO 67 

ABGI.E aZTHEBii aAGd VECTOB ABO MAXH FXBU) • /] 


IF (JORI.Ea. ij 
,, WHITE (6,63) 

53 FOoBAT ('3',' 

00 65 J«1,H0 

I*(J-1)^*3 ♦ 1 

CAU, aAGVBC (J, I, P (X> 1 ) ,P (X) , P (X ,P ax , PI, PO , SIGPaX,SIGaAG) 

HOBTB 

PX»P (I»l) 

EAST 

?T»P(X) 

00 WH 

PZ»PtI*2) 

ia»saSr (P x*px*pi*pi>P 2 *pz) 

l,PX,Pr,P2,ia,PHX,PI,PD,3IGPHX,SiaaAG 


BBIXE16,64] I,PX,PX,PZ,Ra,P 
o4 F0BaAT(d,I3,2Z,9(f7.2,3X)) 
55 CONTiaaE 


57 COHIXMOE 


303 


**io loiaiiV^' Ixio-'TOO auca 
4 1 FOBBAT ('3*,' laPOT BAOXi 


IF (JCOBPH.BQ. 1) CALL COBXPB P,0C,BP,1) 
COIlTIiiaE 
BE TUB a 


DATA RBJECIE0******F1T AEOBTEO') 

BAOXAJ. COBPOBEiil') 

• 3>:i XHPUT SOUTH COaPOBEHT*) 


42 FOBBAX 

43 FOBBAX 
4 4 FOBBAX 

45 FOBBAX 

46 FOBBAX 
SXOP 
EBO 

SUBBOOTIHS SLOCKS , 

coaaoa/SLK/iiflBE(5) ,xa (5,50) ,n (5,50i ,apsa 
GIVES COOROXHAXis OF PBISBATXC BLOCKS IB SB 


*3',' xapui EAST COB poa sax') 

•3' , • isPUT Aaoa xa ip' ) 

<)',• laVESI TO SOUBCE BAGaiiaOES OHLr*) 
>3',' XaVEBI TO SOURCE CUBPOaEaXS' ) 
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: EiLATIVS ro CiJiTEB POIST AI 37'ii,B5m 

3 F08MAT (215) 

. HEAD i12,1) -£SH 

1 FOMAI (i5j 
DO 2 

H£AD (12,5) fl,li 
a-a-1 
iiaB3 (I) a a 
DO 3 j-«l,a 

3 HEAD (12-4) AtAI-ALOH,XX (I,J),2Y(I,J) 

4 FOBBAf (2F6. 1-2F10.2) 

BEA0(12.4) CLAZ,C2.0H,XXA,XYX 

2 coireiaae 

DO 6 lal.apsa 
a>aaBB af 
UEITE Is, 5) i,e 
DO 7 j=i;a 

UCI,J1« (XX (I,J) -34.61) •7.5*25.4 
Il]|, J)* (YY (i, 3) -9,52) *7. 5*25.4 
, iHlTE (6 'ai ixa,J) ,YY (I,J) 

3 FOfiaAT (2f 10.(3) 

7 CONTiaUE 

i CONZiaOE 
HETUBH 
£80 

SUBfiOaZDiS FLO (QU£,ZaO,£)lBB£) 

DATA i«/1 / 

CALL FOG (1,a,0,OaE,Z8Q,i;UBEZ,19l>d.,50,L,A1,A2,A3,A4) 

CH£a-A3/A4 

TS0»-A1/A4 

XUB£E*-*'A2/A4 

L*0 

EEXOBa 

£80 

SOBBOOTIHE FDG (J. BB, 8SZI,OL AX , DL08G , 0,Xd, 8B X,L, Z , Y, A, F) 

**«*««««**«*• «*«*«*«*•« a««« ««•*«« *«««**«•«« ««««««« 

J.EQ.0 £8?UT£ XAXIXOOE i O^ALIIZIIDE (KB) BELAXIVE XO ELLIPSOXO 
Igeooexxc COOEOXBAXEE) 

J.E&.0 DOXFOT FIELD C0BP08B8I5 80BXH ,£ASX, TE8XXCAL 
X8 G£C0£XIC COO BOX a AXES 

I. 8E.0 LAX.&L08G X8 SF8EBXCUL CQ0B0X8AXE5, y*a£OC£8XBXC BAOXOS(KB) 

J. 8E.0 DOIPUT FLO C0BF08BHTS 8UBXE, EAST, /EBXXCAL X8 SPHEaXCAL COOi 

aa.Ey.O use OEFAULX values A£-S37a.16,FLAX*29a.25 

IB. BE. 3 X8PUZ VALUES FOB AE,FLAX OH FXBSX CALL 10 FOG 

HEXI.Ey.O DO HOT BEAD XNPUX VALUES FOB EXIEBHAL FXELO PABABEXEBS 
«B£8 L IS GBEAXEB IHAH 0 

3EZX. EQ.U DO aOX EVALUATE EXXEBHAL FXELO FBCB BUOEL 
HEXI.8B.D HEAD X8P0Z VALUES FOB EXXEBHAL ^ljiLD PABAdBXEES UUE8 
L GAEAXEB 0 

EVALUAXE EXXEBHAL FXELO aOOEL 


HEXX. HE.O 
DLAX 

DLOHG 
i 

HflAX 
iBAXX 
HBAXXX 
HaiXXT 

X. EQ. 0 
X . HE. 0 

tZEBO 

IBAB 


GEDOEXXC LAXIXODE X8 DBGBSES BUE8 J^O 
GE0CE8XBIC LAXXXDOE X8 OECBESS BBE8 


LOHGXXUOE X8 OBGBESS 

liSl 


GSOOEIXC ALXIXUOE 
GEOCBUXBIC BAOXOS 


HBE8 J*0 
HBBH J*1 


BAZIBUB OEGBEE AHD OBOBB OF COHSIAHX lEBBS OF FXELO BOOEL 
" « ■ FiasX OBOBB TIBS ■ n h 

■ ■ « SECOHO " ■ m m H 

" n m IHXaO • ■ a ■ ■ 

FXELO BOOEL COEFFXCXE aXS SCBBIOX HOBBALIZED 
FXELO aCOBL COEFFXCXEHXS GAUSS HOBBALXZEO 

EPOCH XZBE FOB FXELO BOOEL CQEEFXCIEHXS 

BEAH BAOXUS USED XH FXELO BOOEL PUXEMXXAL SZPAHSI08 
(OEFAOLZ « 6371.2) 


aOOBXZ.EQ.a so EXIEBHAL field, solved HXXH BOOEL 
aOOEXX.HE.D EXXEBHAL FIELD SOLVED HXXU BOOEL 
L.EQ. 0 
Ci.GX. 0 
C.LE.0 


EVALUAXE FXELO 
HEAD XH FXELO BOOEL ABO EVALUAXE FXELO 
EVALUAXE FXELO AX OLD XXBE 


BwUXVALEHCE (SHBII (1 , 1) ,XG(1 ,1) ). 
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rAOS 13 

OF POOR QUALITY 


COaflOS /C32?FS/TGi18, 18) 

, 298. 25/ 


r . iSv “*r A 

DATA 1x331/0/ 


DAIA Ai,FLAI/6378.1o,: 

DATA XLA3T/0./ 

DAXA XABAE/6371.2/ 

XP(IFaSI) 110 , 100,110 

JAiXa aADIQS AUO FLAXXaaXHU FACXOa 
USfiD i.ai GtuDEXiC-GEOCEMXaXC CtlOfiJ)ISAXES. 

PAlAaiiias^^^^'^^ XMOEPEHuEax of xhose 

103 IF Wa.aE.Ol BEAD (5, 101 ) AE.FLAI 
131 PgMAX(1X,2F6. 1) ' 


jsaAi i ^rom 11 

ilTEJo, lfi9) AE,FLA 
3 BiJAIj// 3 xt»CcfisTA 
3.3/,22X, 'EABIri dZ 


.BITE 
139 Foaa 
6F8.3. 

iFasx=i 

FLAX=1. -l./FiAX 
E1»0. 

E 2 > 0 . 

E3»0. 

A2»A£**2 

A4«AE**4 

32» (A£*PLAX) **2 

A2B2-A2* (1 .-FLiI**2) 

A434»A4*jf1.-FLiT»*4| 

110 IP - ' 


AE,FLAX 

i, * /*22X, • EwUAXORXAi EABXH BAOXOS 
-.CIP80CAX FLAIXESXaG ',F6.1//) 


IF (1) 19,1.2 
II {xh-riiaf) 17,19,17 

BE Afi J5, 3^ i “?Sf » » aaiXTT .BdllXT, aOOEXT,K , XZBBa, AoAa , 




103 

104 
105 



'. 103 l 

! (20A4 


IF 
BE 

Foaalx 

.FOBEAX (Sx, 'FIEXD dUOBL OEdeb rl 
I,'£XZEBilA 


•IT, SHXTtT-aODEXX,K,XZEao,AaAB 

SST-'ETELn ’inriPT .42.' lu Bir\riic< ■ 


BAOXOS ',F6. 1//) 


SSZ.'FIEXO EOOEX BEAU 

aiXM-O 
X£ap*o. 

FOBflA^'’'( 2 l 3 " 5 Fll* 4 f“"“'®^““'^*“'‘'®^’'““'“^^““ 

IF (B.xa.O) ' goxo7 ' 

aA^-(aAX 0 (S,aAXiij) 


sad 

>aXBB 


iijjjli, 

GXTi«,il 

dP>id 

CB* Am^ < j 

Gla-i,«j sdiiB 


-- JlS.Sl -aXTHB 
i|B^-iaixi <xEai.Afls(Gxaa)) 

IF (B.EQ. 1) GOI05 

G(»ii.ar»ali“ 


7 

106 


107 

102 

3 


GO ro 5 


Bxxaa 



Gift {a.BiiGiiTia 
IFja.Ew.1|^ GO_ig 106 


Gift (fl-1,i() 
GO xd 10& 


•axTiaa 


COHTiaUE 
IF ■ 

FOi 




S1,E2,E3 


00 12 B«2,BAXa 
DO 12 d*1,a 
fli»a - 1 

IF (B.E3.1) GOXO10 

iFjil.ax.aaixii) ■BxiE(6,9)a,a,G ta,a) ,g iax,a) , 


59 


nnnnoriD 


BisiyEss AMD Techmological Systems, Inc. 


10 

11 

12 

103 

13 

14 


15 


16 

17 


210 

220 

230 

273 

240 

250 

260 

13 

19 


C 

20 


OOT (3,a),ar (3,9) ,0TT(9X,ei) 

1 ? (3. Xi:. 61X111] ii&XT£ (6,9) 3, 9,0 (3,9) ,0 iHI.M) , 

^FO^liT*5l3 V ,SII(9I,3) ,OXTT (3,9) ,GITT19I,3) 

Go“lO ^ 

CO BIX 3 OE 

IF (3.GT.39XTXT) ■filT£(6,11) 3,9,G(3,9) ,GI(3,9) , 

I? (B.ll. H9XXTI] UBXI£(6,11}3,9, G(a,9),GZ(3,3) , 

SGITJS.d) .GXTT (i.ar 

FOMAf (iX3,F1 T-4,1 1i,F1 1.4,1 U, F1 1.4,1 1X,F 11 .4) 

COHXXBUS 


IF (900EXX. 3E.01K6XXS(6«103) El, E2,£3 
F0a96X (//SX.aEIXXFLO, ,3F10.2) 

FOB9XX ilBl) 

IF (XEBP.Eu.O.) L*-1 
IF IS.SE.O) GOI017 


SH9II (1,1) =-1 
OO 15 3> 2,94X3 


(3,1) »SH 
(1*3) >0. 


SH9IX (3-1,1) ♦FL01I(2*3-3) /FLO AX (8-1) 


S39IX (3,1 
SH9XI “ 

JJ>2 
DO 15 9« 2 . 3 

3H9II (3, af »SnaiX (3 ,9- 1 ) ‘Soax (FLOAT ( tB-9» 1 ) • J J) /FLOAT (3 ♦ 9-2) ) 
SH9XX (9-1,3) >S39XZ (8,9) 

Jj3l 


DO 16 3-2,94X3 
00 16 9-1,3 

G ( 3, 9 ) -G (3 , 9 ) « SH9XX (3 , 9) 
GljB,9) _»GI (B,a)-SH9I 

IF . 

IF (9.EQ.1) GOT0 16 


(3,9) -G (H,B) -inaxr (3, 9) 
r (3,9)«GX (3,9)«SH9XX (3,9) 
ri(3.91 -GXT (B,9)*SHaXI (3,9) 

F (BalTTX.GI.a.A30.3.L£.d) GITI (3,9)-G 


GXTX (3,9) «SH9XX(3, 9) 


G (9- 1 ,3) -3 (9-1,3) •SH9IT (9-1 , 3) 

GT Id- 1 ,3] -(JI (9-1 , 3) •SaaXX (9- 1,3). 
aXX(a-l,B) -GIT (9-1,3) •399X1(9-1 ,3) 

IF (BaXXfx.GI.0.A30.8.LE.3)GXTT (a-1,3)»GXXT (9-1,3) •SH9IX (9-1,3) 
COBTIBOE 


x-xa-xxsao 

00 18 3-1,94X3 
00 13 9-1,3 
XGX-0. 

XHX-0. 

I?(9.Ba. 1) GO TO 270 
IF (3. GX. 39 XX XT [ GO XO 210 
TGX-GTTI (3,9) ‘i 
XHX-GXTI (9-1 ,B)*T 
IF (B.GT.3B4XlZ) GO 10 220 
TCX-(XGX ♦ GXIJB,ai)*X 
THX- (THX ♦ GTIia-l.ar)*! 

IF (3.GX. 394XZ) GO XO 230 
XGi-(XGX ♦ GX(I,9))*T 
XHX- IXHX+GT ( 9- 1, 3) ) *t 
XG3-XGX>G(3,9) 
iai-THX*G]a-l,3) 

XG (3,9)-XGX 
XG (9-1,3) -THX 
GO XO 18 
CO BIX 3 OE 

IF (3.GX.39XITT] GO tO 240 
TGX-GXTT (3,9)*1 
IF (3.GX. BalXXI) GO TO 250 
XGX- (XGX^GXT ( 3 , 9 ) ) *T 
IF (3.GX.394XZ) GO XO 260 
XGX- (XGI4'GX(B,B))«T 
XGX- XGX >G (3,9) 

XG (3, 9) -XGX 

C03XX3UE 

TLASI-TH 

DLAIH-DL4X/57. 2957 795 00 
3I3L4«SI3(0L4XlI 
BLOHG- DL03 G/ 57 .2 95779 3 00 
CP9-C0S(BL03G) 

SPH-SX3(aL0BGi 
IF (J.EQ.O) G01020 

Q IS GSOCEITBXC B40XOS 3 9E3 J-1 

a- a 

CT-SI3X4 
GO XO 21 
SI3L42-SX3L4«*2 
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z^ 

22 

23 


9 

10 
3 


u IS s^oi£ixc zhsu j=o 

iiT=>i 

C3Si.A2-1.-SXIiLi2 
DZ32*i2- A2B2«iiliLA2 
DJia=Si5T (D3N2) 

,<0SH)+42)/({u*0EH)*B2) 1 
i/S UE T (IAC«COSXA2 *SIE I.A2 ) 


Poor 


iO 


Quality 


42)/({W*0EH)*B2) 1 **2 
A2*SIEi.A. 

•0£E) * {&4-&4B4«SXHi.A2) /0EE2) 


lOMS 


F4C=»ll(^* 

CI»3Xi*l/ 

B»saET(a*iw«-2.'' _ 

3T*SwaT11.-CX**21 
N3 iX^aiS OilUli, SAaE) 

NEXIFxKEXI 
CALX. EA3F 
la BP 
F*a 

XF |J) 22,23,22 

2a *o R 
fiXZUBE 

lEAESFOBflS FIXIO TO GEQ03XXC OXBECX 
SIEO=3IELA«SX-SwEX (COSLA2) «CZ 
COSD*3wEr(1.0-£IMD**2) 
ja-fli*c3 S3-oB* SX HD 
Z=aT*SXHfl-ai*CCSD 
SEXaEH 
END 

SOBEOaXXHE HAGX 
COtiaCH /C32FFS/G(13,18) 

COdaOH/FLDCO a/SX.CX,SPH,CP3,B,2HAZ,aX,BP,BE. a.ABAE.EI, £2, S3.HSXX 
OXaZHSXOH ?( 18, Id) .DP (13,13) ,C0H3Z (13, 13) ,SP ( 18) ,CP(13) ,FH (18) ,FJ ( 
1 18) 

XF (P (1,1) .Sw. 1.0) GO TO 3 
P(1,l)»1. 

DP (1,11-0. 

3P (11=0. 

CP(lj»1. 

DO 2 H=2,18 
FH (H)=H 
DO 2 8=1, H 
FH (81=8-1 

CO ESI (H,8) -FLOAT ((H-2) **2- (8-1) **2) /FLOAT ( (2«H-3) ♦(2*H-5) ) 

SP (2)=SPH 
CP (2) =CP8 
DO 4 a=3,H8AX 
3? (31 =SP (2 ■ 

CP ]a) =cp i 2 

AOB-ABAE/B 
AE=AOfi««2 
BT=0. 

BP«0. 

3E=0. 

DO 3 3=2,E8AX 
AB-AOB-AE 
DO 3 8=1, H 

? 

0 


j *CP 

(8-11 ♦CP(2j 

•SP 

(8-1) 

) *CP 

Ih-i 5 -sp {25 

♦SP 

(8-1) 


0 3 8=1, H 
F (N-ai 6,5,6 
(H.H) =ST*P IH-1,H-1) 

P (HgH^=3r*DP (H-1.H-1 


GO , 
P(H,8) 


1) *CT*P(H-1 ,H-1) 

CI*P (H-1,8) -COiSI (*,3)*? (H-2,8) 

HOXB : CCHSZ(2,1)=0 

SX«P (i- 1, 8) -COH3I (H , 8} « DB (H- 2 , 8) 


DP (H,8) = cr*DP(ff-1,8) 

PAB=P (11,8) «AB 
IF (8.20.1) GO TO 9 
X2aP=G (H ,8) *02 J8) -»G (8-1,8) *5? (81 
3P=aP-(G (H,8)*SP(8) -G(8rl,H) *CP (8).) 
GO TO 10 

X2aP=G(H,8) «CPJ8' 

»8T*TESp*DP 


•FH (3) ‘PAH 


aT=8T*TEap5»DP (S,8l *48 
3B=8B-TE8P*FH (H) «PAB 
ap=BP/sr 

XF (HEXT. GT.O) CALL EXTFLD 
B=S0BT (3T*3T>BP*BP>BB«BB) 

BETOBH 

EHO 

SOBBOaiXHE DATA 

COaaOH/DAI/B (60,31),X (e0,31) ,P (60 ,3 1 ) ,B (bO ,3 1 ) 
DO 1 X»1,31 

- . 


aSAD ^14 12) (3(J.X),J=1, 
BEAD (16,12) (T(J,X),J-1, 


DO 3 


■1 ,31 


aO) 
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% 

12 


B2A0 (13, 12) i? (J,I) ,J:il,bO) 
DO 4 1=1,31 

3240 (20,12) (E (J,I) ,J=1,oO) 
FOEdAI (1127.2) 


BSTOBH 
EBD 

Sa3B0UIX8£ 
OISE HSiON 

coano» ■ 

COflaO»/BBC/P 
•0FDP4(14) 
B24L*S D, 
aBlI.*4 _ 
iI2iC«70. 
D£C=-3. 

PB»0. 

PT=0. 

rp=4. 

YC»0. 
DZ=Z2-Z1 
i1f»Z1 - z? 

,,2?= ^2 " mF 

DO 57 n»1, 
MBIT2 (6.40) 


JUNJJOPTl 

5w),2(50) ,X(50) 
bil 

OFDP3 (14)., 

D.uFO$1,OfOP2,OFOP3,&FOP4,OFDP,03,DP,S3IlD,OC,i)CC 
L2, 3 , B 



43 FO£aAj(2l3,2F3.2) 


BD 

XX,BPSa,XF,XF 


L»(II-1)*3 


HBP=lia£fi (IZ) 
1IB1*HBP«’1 
II 
XX 
DO 

ai (t)=iu 

■II 


,XI-»B1)=TI| 
SiBP 





Siiiu 


BP 
XX 


XX, 1) 
XX, 1) 


C WiXXE (6,50) ,1-BBP, 
S3 FOfiaiZ (2X3,5X10. If 


X*1) -IX(XI,X1 
X>1)-II(XX,I| 
x (iJ*Bxlx) ♦ai (X) • 
x,ii6p,xi(ix,xi-ij , 


•BY 

XX 


{iU 


x),as (X) ,xi(xi,x*i) ,ii (XX, X) 


COBTXHUE 


/BS|X| 


/BS| 


DO 2 x=i,aai 

1(1) »XX(IX,X)-XF 

x]xL»ii(xx,x)-if 

CONTiaOE 

71=0. 

72«0. 

73=0. 

74=0. 

75=0. 

76=0. 

DO 3 X=1 ,8BP 
X1=X X) 

X2=X X«1) 

n»i X) 

12=1 X*1) 

ox»Bxr- 


0I> 

05^ 


31 

■as 


al=SCBT(X1* 

a2=s3iTjx2* 

ai 1=3021 (B1 

a H a— Tr* fw 7“ « 


xi*zi*n*ii) 

“'•X2+ 12*12) 

_ soar (B1*B1*Z1F*Z1P) 
ai2=SQBZ i1«E1'»Z2F«Z2F 
£21=s3bX B2«B2*Z17«Z1F 
B22»SwBT ;B2*&2*Z2F«Z2F 


01= (x't*DX*X 1*01 
02>(Z2«0Z*I2«0I 
PP*iX1*I2-X2*ri 

Ti i»aii*oi 

T12=B12*01 

T2 1=B2U02 

T22=B22*D2 

71 1=ai1*ZlP 

F12-B12*Z2F 

F21=B21*Z1F 

722=B22*Z2? 


/OS 

/OS 

/as 


47 



£.■ 

FOaaif (12ir8.2),2I3) 
7»AiOG(t22*7 11>(712* 


712,721,I12,I21,Z1,Z2,PP,BJ1,ii22,i12,B21,X,aap 


. , F21 

g^iiOG (I22*T11/(T12*t21 

X7(?P.22.0>) SC ZO 11 


i! 
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a=»rAS i22F*D2/ (PP»322) )-AIiN (Z2F*DV (rP*ai 2) ) - 
♦ATAi* Uli;*D2/ kP 5*321) ) ♦AXAS (2 1P*D1/ (PP*H1 1) ) 

1 1 CO tiTXii JS 


12 


31 


707 

57 


71=71*SC*F-C2*li 
72»V2*SC*li+C2*F 
73»73*C(I) *0 
74*74-3C*F-S2*H 

75- 75-3(1) *t 

76- 76*n 
CO STIHUS 

IF (J0PT.3«. 1) SO TO 12 

FOa ZSTLAkTlsa SOOaCE COaPOMEMTS 


)d 




-75 
-73 
-76 
-72 
= 71 
-73 
=74 
= 72 
= 75 

=LL*72*a*74*«*75 
-Ll*71*a*72*S*73 
, =U.*73*H*75*S*76 
•73+P (L*l) *75*P (X,*3)*76 FH 
•72*P (L*1)*74*g(L*'3)*75 ♦ FP 
•71>P (L*1)*72+P(i.*3)*73 ♦ FT 

•0FE54 (L*1) (i*2) *0F0P4 (1*2) *P (L*3) *DFDP4 (L*3)*IC 


DFDP1 (i*1 
DFDP1 i*2| 

DFDP1 1*3 
0F0P2 L*1 
0F0P2 1*2 
DFDP2 1*3 
DFDP3(L*1 
DF0P3 i.*2 
0F0P3 L*3 
0FDP4 l.*1 
3F0P4 L*2 
0F0P4 L*3 
F£-P (L*2 
FP-P 1*2 
FT-P L*2 
IC-? L*1 
GO ro 31 
coraiauE 

FOB ESTIBAXXHG SOOBCS aASKITOOBS 

(ii*73*B*75*8*T6) 

’lL*72*a*74*M*75) 
i.L*71*a*72*8l*73) 

- - lL*71*H*72*a*7 3) *a* (LL*72*a*74 *N*75) ♦ 

♦N*76) 

*OFOP1 III) *FB 
•0FDP2 |llj *FT 
•DFDP3 II) *F? 

•0FDP4 III) *IC- 

2aifS'(6,707)_ U-OFDPI ill) ,OFDP 2 (II) ,aF0P3 (II) ,aFap4(II) , 
U..at2,71,V2,75.74,75,76;F,a. »‘PP 
FOBaAI(I4,4F7!2,5F6r2,6F6;2i4f5ll) 

COailHOE 
BEIUfiH 
EHTBI AEGI 
AiaC-70. 

DEC— 5. 

ABC-. 0174333 

U.-COS (AiaC*ABC) «COSiOBC«AHC) 
a-CJS (AIBC«ABC)«Sla (0£C«A£C) 
a-SZV(AIHC«A£C) 

BEZOBtl 

EHO 

SaSBOOZIBE HAG7EC (J, X,PX,PXr PZ, BttI,PX,PD,SISPHZ) 

BEA1*8 OHOiaxni 
DiaEBSIOa EC07 (3,3),A(3) 

coaaoH/oaAi/oaoBax 

COaBON/POS/aO,XLAX(4) ,ZLOa (4) 

^^^LOC|li^JJ, aoiil) » (JJ-1) *8018 - (JJ**2-JJ) /2 

FI-70. ^ 

FD— 5. 

P»SQ£Ii?X**2*PI**2*PZ**2) 

Pi-lBsia (PZ/P) 

PO-ABSia ]PT/ (P«COS (PI) ) ) 

FZ-SIH ((FI)* (ABC)) 

rT-COS{(FI}« (AEC)] *Sia i(FO)« (ABC)) 

FX-SQBTi 1. -F X**2-FZ**2) 

Pai»lBCOS( (FX*FX*fX*PI*PZ*PZ)/P) 

PI- (PI) / (ABC) 



* II 


PO-(PO)/ (ABC) 

B2-PI/P 
B3-PZ/P 
CBE-XIAI (J) 
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995 

23 


993 

995 


U 

12 


15 

903 



I20=YiOa(J) 

rfl£BE*6371.2 

XP=I 

IP3=I ♦ 2 

DO 20 IC=IE,X?3 

DO 20 JC*If;i23 

LC^LOC (XC,0C,S0) 

iF{ic.i2-jq i.i;=xoc (Jc.ic.sD) 

uaxrx(6.99£) ic, JC.Lc.DaoadXiLQ 

yoEaArj* bcov Hiiaix *- 3 x 10 , 020 . 12 ) 

sail 
HaiX2 

FOaaAlC* 'aiovxc .^90’ '•,3X10,2X8.2) 
FOaaAT pox, 8215.3) 

CAXL FL07ck£,I20,IHa22) 

COSX® (-Oa2*P2 - THO^PX ♦ TURSE^PY) /P 

sxa?*scaili. - cosp*cosp) 
i(1)*(fHi2I - COSP*B1) /l?«SiaP) 

-IBO - COSP*B2)/{i*SXHP) 

•0a2 - C0SP*E3j/lP*SXHPi 

* 0 . 

DO 12 X»1,3 
00 14 iC-1,3 

axsrai^siopux * i(K) «Pcov (a,D «A(L) 
co«iiao2 

SXOPHX=SuBI (SXOPHX) /. 0 17 4533 
* '-•81 
•82 

, - 

CoaxG^o. 

DO 16 X«1,3 
DO 16 K«i;3 




siaaio^siGaAa ♦ aiki •pco» (K,l) *4 (D 
axr2ib-900i J,x.cdsP,sxopili,slGaiG 

omii (10Z.* aAG727.900 *,2x10,3X10. 


aaxr 

POkRk^ 

B2XOaH 

2HD 


X) 


SOBBOaTIHE COBI,PB(0,S,80a,8Z) 
a£Ai«8 D(l) ,S{1) -tlovaiH 

s^>pza xa 0 x«ttoa pta j>cox 


330 

323 


_ •COi PTE 

0 AaaAX HOLDS HOaHAL EQUATIOMS OE CUTAEIAHCj 

s xs pa za tool aebax 

12(81.20.0) GO TO 320 
DO 300 iiSfi 
LC-iJ-1) *aoa-(j*j-3*j)/2-i 
DO 300 a » o,Hoa 

K »1 


aaiTS ( 6 , 72 a) x 
U£XT2i6.73Ul 0(8} 
DO 40d J«2,aoi 


35 0 K 


K>J 

DO 350 X>1,J 
S il]«D(S) 

”- KVHOa-t 


aaxZE (6,720) J 

, aaiT£ (6,730) (SII),Z»1,J) 

43 0 coaiiaoE 
BExaaa 

720 POEHATt/* 0 »*,X5) 

73 0 FOaHATi* *,13210.2) 

Eao 

SOBBOUTZaE ZSXlT(LZ,aa,A,i) 
DOUBLE PBECXSZUa 0PXV,D 
ZOXGL>0 
LTBOH«1 

Z2(LL.LZ. 1) GO TO 900 
LL1*ti-l 
R 1>0 

La>aa-LL 
iao>-La 
DO 90 a>1,LL 
IHD « ZaD>Li 
xPZT ^ zao^i 
L £ H 0 > R -1 
ZOL>A (KPZVt 
DO 30 I>X,Ll 
iao>zaoEi 


sua,A2,aii).,A(i) 


SAXazz 
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10 


20 

30 


UO 


> 0 


70 

30 

30 


110 


120 

130 


0SUd>0.O0 

30,30, 10 

LIMO-I-K 

DO 20 L>1,L£il0 

osua>osaa«’A iLiii;}«i(LA!i?*LXNO) 

LAil?>LAli?«'aB->L 
CO NTZIiOS 

DS0d>A (XHD)-0S0{1 
IF (I. aZ. SI GO IC 70 
IF (DSOa) 900,900,40 
CONTI HOE 

lOIG^ALOGIO (TOX/SIIGL (DSUai ) ->.5 
IFTIOIO. I.E.iOIGL)aO TO 60 
XOIGl.>IOXG 
XISOB>X 

DFXV>0S0Sr(DS01i) 

A1»(1.00/0PIV) 

A2» (1. 00*Dai.E (A1) «0PXT) /DPI? 
A(IllD)>OPIT 
BiK) »I>?IT 
GO TO 80 

A(xiiO)^A2*osaa«oBX2(Ai) «osua 

CONTINUE 
CONTINUE 
DO 152 S>1,LX 
DPIV-A iKPIVI 
A1»(1.D0 /DPIV) 

A2*M-D0-DatE(A1) •0PI?)/DPIV 

AfK?I7)»A2*0E££(il) 

i]lL-K+1)»A(KPIV) 

LEND«5-1 

IF (LEND) 130, 13C, 110 
DO 120 l»1,t£NC 
IND>KPXV*^L 

C^NIlioi ♦OfllJ (A1) • A(X« J) J 


OF POOR QUAU1Y 


152 


140 


150 

151 

152 


LABFlsA2*DSgB+D3LE(A1) •DSOB 
(LENDpSl, 1*1, 140 
150 L»1,l£aC 


150 

170 

130 
921 
93 0 

92) 


IF (S.EU. IL) GC TO 
XND»KPI7 

K?IT«KPIT-XH-1-K 
IANF»XNO 
DO 151 X«K,I.X1 
LANP«XAJIP-La-I 
DSUa>A(i.ANF) 

AaABF)sA2* 

IF (I.EN£ 

DO 150 i.-., 

iino>i.anf«L 

S^UNO|j;A(XXNO)»DSUa*A(XNO*X) 

CONTINUE 
CONTINUE 
DO 180 K>1,LX 
LIN D»KP 17-1 
LABF-KPIV 
DO 170 X>K.LL 
DSUB»0.D0 
00 160 l.-KPI7,IN0 
LI ND>LXND*'1 

DSUB«OSUB«’A (L) *A (LIND) 

CONTINUE 

AIKPX7)»DSUa 

LX NO-LIN D«’ LB 

KPX7-KPX7* 1 

CONTIBUE 

B(K)»0.D0 

5PI7»KPI7»La 

IND-IND^BB-K 

CONTINUE 

BBXTEi6,921) XDIGL.LTBOW 
F0EaAT(//2X,' •*•••• ISIN7 ••••••• ,2X, 2X7//) 

REXUBH 
XDXGL— 1 
LTBOB-X 
aBITE(b,9201 
FOSBAIiSl, '* 

STOP 13 
BETUEN 
END 


LTBOH 

• « « 


IN7EBSXON FAILED AT BOB', 16) 


/• 

// eZEC OLINXGO8,E£GION.GO-500S 

//GO. FT05F001 DD DSH-F9#BG. GaCOEF(PO3Q0272),0ISP-SilH, LABEL- (,, , 
// DO DSN- ICOaa.INV2BT. AfiEA(LqC60Z32) ,0XSP-SUB, LABEL- (, 


IN 

9 9 
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0SJI»]fIHAil.U3-4££A-ao.3ATi,DXSP»SHB,LAB2i3 (, , ,XS) 
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Table 1 


Table 


Table 3 


TABLE CAPTIONS 


Summary of regional prismatic model computations using a simple prism 
model (Figure 10) for the Kentucky body geometry. 


Summary of regional prismatic model computations using a detailed 
prism Figure 11) model for the Kentucky body geometry. 


Summary of Mosaic Dipole array model computations using the regional 
geometry shown in Figure 17. 


67 



Bisiivess a.\d Technological Systems, Inc 



792.0 


750.1 


comvo csjmco 

O ^ • rn ro 


CM ro 


O cn 
CO ^ 


o 

d 


^ CM 


O 


— • O 


O «-• CO 


CM O 


O CO iT) 
^ CO 


O uO 
O VO 


CJV 

O 

VO 

CM 

LT> 

50 

CM 

r«». 

fs. 

en 

<.X5 

CO 

VO 

CM 


CM 

CO 

ir> 


VO 

CO 



LT) 

CO 


m 

03 


^csjcn ^CMCO ^CMco ^CMro 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


68 











BisiyEss AMD Techmological Systems, Imc. 



69 


334.3 










Blsimess and Technological Systems, Inc 




BuSiyESS AMD TECHNOUHilCiL SYSTEMS, INC 


FIGURE CAPTIONS 


Figure 1 Equivalent source representation of the magnetic anomaly field at 
height of 325 km derived from Magsat data. Units are nT. Albers 
equal area projection. 


Figure 2 Apparent magnetization contrast in 40km thick layer, obtained from 
Magsat data. Contour interval is 0.1 A/m. Albers equal area 
projection. 


Figure 3a Bouguer gravity in the vicinity of the Kentucky body. Contour 

interval 6 mgal . Refraction profiles in fence diagram form from 
Warren (1968); depth scale marked off in 10km intervals. Inferred 
position of Grenville Front in heavy dashed line. Light dashed 
line is aeromagnetic low from Figure 3b. From Mayhew et al (1982). 


Figure 3b Aeromagnetic anomaly contours in same area as Figure 3a. Values 

are hundreds of nT, contour interval 400 nT. Relative to arbitrary 
datum. -30 mgal contour from Figure 3a shown. From Mayhew et al 
(1982). 


Figure 4 Geometry of model Kentucky body; three regions discussed in text are 
indicated. Angled boundary shown indicates area of Figure 3. 


Figure 5 Computed Bouguer gravity due to model Kentucky body. Contour 
interval 6 mgal. Compare with Figure 3a. 


Figure 6 Computed magnetic anomaly due to model Kentucky body with arbitrary 
datum shift. Contour interval 40(1 nT. Contour values are hundreds 
of nT. Compare with Figure 3b. 


Figure 7 Magnetic anomaly due to Kentucky body at satellite altitude. 
Contour interval 1 nT. 
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Figure 8 Light broken line is one contour line selcted from aeromagnetic map 
of Zietz (1982) to indicate extent of highly magnetic source 
region. Heavy solid line is -30 mgal gravity contour to indicate 
extend of Kentucky body (KYB). 


Figure 9 Tectonic elements in region surrounding Kentucky body. -20 and -30 
mgal gravity contours from DOD compilation shown as heavy solid 
lines to indicate significant highs (h). Cincinnati Arch delineated 
by zero level structure contour (dot-dash line) on top of Trenton 
(uses and AAPG, 1962). Heavy long-dash line is inferred position of 
Grenville Front. Generalized faults of 38th Parallel Lineament 
principally from U.S. Geologic Map and Ammerman and Keller (1979). 
RCG is "Rough Creek Graen" (Soderberg and Keller, 1981). Other 
symbols are as follows. WL = "Woodward's Line", ECGH = East 
Continent Gravity High, MMGH = Mid-Michigan Gravity High, LFZ = 
Lexington Fault Zone, JD = Jessamine Dome, PMT = Pine Mountain 
Thrust. Position of Kentucky body as delineated by gravity contours 
labeld KYB. Aeromagnetic lineament shown as short dash line. Small 
circles are selected basement core locations. Solid circles are 
medium- to high-grade metamorphics. Open circles are felsic 
volcanics; circles with dots are basalts. Core samples of low-grade 
metamorphics, sedimentary rocks, and plutonic rocks not shown. From 
Mayhew et al (1982). 


Figures 10-12 Blackened areas refer to first, second, and third model source 
regions, respectively, referred to in text. 


Figure 13 Magnetic anomaly in the total field due to source region two 
(Figure 11) computed at 325 km for unconstrained magnetization 
di recti on. 


Figure 14 As Figure 13 for magnetization direction constrained to be in main 
field direction. 


Figure 15 As Figure 14, but for region three (Figure 12). 
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Figure 16 nifference between data of Figure 1 and data of Figure 15. Map 

shows anomaly in the total field without the effect of the extended 
source region. 


Figure 17 Geometry showing dipole locations and mosaic regions I (•)» H (+) 
and III (a) for the dipole array models. 
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FIGURE 4 


77 




FIGURE 5 


78 





FIGURE 6 
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FIGURE 7 
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FIGURE 10 
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FIGURE 12 
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FIGURE 15 
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